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Abstract 

We compute the quasi-potential and determine minimizing paths 
for an action functional related to scalar conservation laws on an in- 
terval with boundary conditions in the sense of Bardos et al. ([3]). 
Taking as input an exclusion-like flux function /(p), a strictly convex 
entropy h(p), and boundary data pi,p r , we obtain a generalization of 
the functional derived in [23J for the stationary large deviations of the 
asymmetric exclusion process. 
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1 Introduction 

Let be given K > and a flux function / on [0, K] satisfying 

feC 2 ([0,K}), /(0) = f(K) = 0, f"<c<0 (1) 
We consider the scalar conservation law 

d t p(t,x) + d x f(p(t,x)) = (2) 
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on (0, 1), with boundary data 



p{t,0 + ) = p h pit, I ) = p r 



(3) 



where pi, p r G [0, K]. It is well-known (see e.g. jB]) that, for given Cauchy 
data, ©) has in general no classical solution, as even for smooth data, shocks 
appear in finite time. On the other hand, uniqueness fails for weak solutions. 
For the problem on R, the unique physical solution (called the entropy so- 
lution) is selected among weak solutions by the so-called entropy condition. 
For the initial-boundary problem (J2J) (J3]) , there is in general no entropy solu- 
tion satisfying (EJ) in usual sense, but a unique entropy solution that satisfies 
03]) in the weaker sense of Bardos et al. ([3]), hereafter called BLN boundary 
conditions. In short, pi and p r in fl3]) must be viewed as sets rather than sin- 
gle values, and the definition of these sets depends on the boundary location 
through the outer normal. This phenomenon arises from the formation of 
boundary layers in the inviscid limit for the viscous parabolic approximation 



In statistical mechanics, the asymmetric exclusion process (ASEP) and re- 
lated interacting particle systems in the sense of [35] are microscopic models 
with a single conservation law, whose large-scale behavior (the so-called hy- 
drodynamic limit, see [30]) is governed by © . For many such models, see e.g. 
[13] , it is established that the empirical particle density converges in probabil- 
ity to the unique entropy solution of ©. It was also established (pfl[2]) that 
boundary conditions of [3] are obtained if the particle system is coupled to 
reservoirs at the boundary of some open domain. Let us consider ©-© on 
(— oo, 0)t x (0, l)a.. We introduce an action functional I[p(-, ■)] that quantifies 
deviation from p(., .) being a BLN entropy solution. This functional is be- 
lieved to be a dynamical large deviation (LD) rate function for the time evo- 
lution of the empirical particle measure in ASEP-like models. Our goal is to 
compute the quasi-potential V[p(.)] for a function x G (0, 1) p{x) G [0, K], 
defined as the minimum energy of a path leading from a stationary state to 



and determine the minimizer(s) in ©, i.e. paths followed by the system to 
create a fluctuation from expected stationary states. In ©, S denotes the 
set of stationary solutions of ©-©. This set contains a unique constant 
function p s (.) outside the phase transition line pi < p r , f(pi) = f(p r ), where 



of©. 



P(-): 



V[(p(.)} := mi{I\p(., .)] : p(-oo, .) G S, p(0, .) = p(.)} 



(4) 



2 



it contains arbitrarily located shocks connecting p t and p r ([2]). An analogous 
picture was established in [2T] for the stationary state of ASEP. This picture 
was extended heuristically to more general models and truces in [12], and a 
rigorous proof of this extension given in [2] using BLN boundary conditions. 
Under suitable assumptions (see e.g. in the context of lattice gases), if 
one can prove the dynamical LD principle with rate function /[/?(., .)] for the 
empirical measure under the stochastic dynamics, this implies the static LD 
principle 5'[p(.)] for the empirical measure under the invariant measure of the 
stochastic process. 

The quasi-potential was introduced in [25] in the context of stochastic per- 
turbations of dynamical systems, like for instance 

dx\ = -VU(x £ t )dt + e 1/2 dW t (5) 

where W t is a standard Brownian motion and U (x) a potential on R with a 
single minimum at (say) x — 0. In this case, instead of (J2]) - ©, the limiting 
evolution problem is the ODE 

x t = -VU{x t ) , (6) 
the dynamical action functional for a path x, = (x t , t £ (— oo, 0]) given by 

= \ f [x t + VU(x t )] 2 dt, (7) 

and the quasi-potential V(x), where x £ R, reads 

V(x) = inf{I(x ) : x_oo = 0, x = x} (8) 

A standard computation leads to V(x) = 2U(x), in agreement with the here 
explicit invariant (and reversible) measure // = Z(e)~ 1 e~ 2e u ( x 'dx for (jSJ), 
which identifies 2U as the stationary LD rate function. The same computa- 
tion shows that the unique minimizing path (or "fluctuation path") x, in (jSJ) 
is obtained as the time reversal of the solution of ((6]) starting from x. This 
is a particular case of the Onsager-Machlup principle, which states that for 
reversible dynamics the fluctuation path is the time-reversed relaxation path. 

Our problem is motivated by the study of nonequilibirum stationary states 
(NESS) in exclusion-like models coupled to reservoirs. A general feature 
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of NESS is the occurence of long-range correlations, one signature of which 
is the non-locality of the stationary LD rate function. The first example 
of such a rate function was obtained in [22] for the symmetric exclusion 
process coupled to reservoirs. The same authors ([23]) later considered the 
asymmetric exclusion process, which shares the non-locality feature, but is 
somewhat more complex due to the presence of two distinct regimes p\ < p r 
and pi > p r . The weakly asymmetric exclusion process was treated in the 
same spirit in [20] in the regime pi > p r . The result of |23j is our primary 
motivation. In [6] was introduced the so-called "macroscopic fluctuation the- 
ory" . One outcome of this theory was to recover the static functional of [22] 
by computing a quasi-potential and show that the fluctuation paths obey 
a nonequilibrium generalization of the Onsager-Machlup principle, namely 
they are time-reversals of relaxation paths of some adjoint evolution prob- 
lem. Only in the reversible case does the adjoint problem coincide with the 
original problem. This approach does not depend directly on microscopic 
details of the particle dynamics, but rather on the dynamic functional. Thus 
it is applicable to certain other models. It was recently ( [HI E] ) extended to 
the weakly asymmetric exclusion process, whose fluctuations in the vanishing 
noise limit e — > are described formally by the infinite-dimensional diffusion 

d t p + d x f(p) = ud x [D(p)d x p]~e^ 2 d x [a(p) 1/2 B(t,x)] (9) 
p(t,0 + ) = Pl (10) 

P (t,l~) = Pr (11) 

where p = p(t,x), B(t,x) is a space-time white noise, and v a viscosity 
coefficient, with the particular choice of flux /, diffusion coefficient D and 
mobility a given by 

f(p) = p(l - p), D{p) = 1, a(p) = p(l - p), (12) 

The action functional I v for (Q is an infinite-dimensional version of ([7]) that 
forbids violation of ( 1TU|) - (ITT|) . Its quasi-potential has an explicit local expres- 
sion on the torus ([?]) that is valid for the most general /, D, p: 

S[p(.)]= [ h[p(x),p ]dx (13) 
Jo 

where h(p, p ) is the relative version (see ff25|) ) of the entropy h given by 
Einstein's relation 

h"{p) = D(p)/a(p) (14) 
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in agreement with the uncorrelated stady state of WASEP on the torus. In 
the nonequilibrium case (IT0l - (ITTl) with pi ^ p r , for (JT2|) . the non-local func- 
tional of [20] is recovered in the regime pi > p r ([§]) and another one ([9]) 
obtained in the regime pi < p r . 

Our purpose is threefold: to generalize the static LD rate function of [23J, 
identify it as the quasi-potential for the action functional associated with 
©"©J and determine the adjoint evolution problem that yields the fluc- 
tuation path. The variational problem induced here is very different from 
(]9|) . (TTOj) ( TTTj) because the dynamic action functional is no longer an infinite- 
dimensional version of ©. In particular, Hamilton- Jacobi techniques from [B] 
are not effective here. We take as basic input for our problem the exclusion- 
like (that is, satisfying ©) flux / and equilibrium entropy function h(p) 
defined for [0, if]. For the ASEP case studied in [23J, 

K = l, f(p) = p(l-p), %)=plogp+(l-p)log(l-p) (15) 

The action functional is decomposed into a bulk and a boundary term. The 
former is determined only by the pair (h, f), the latter also by the boundary 
data pi, p r - 

The bulk term is based on [281 HHJ S EZ] where two closely related func- 
tionals were introduced as LD rate functions on the torus for the ASEP and 
for © in the limit v = e — » 0. These functionals are supported on weak 
solutions of © and measure the positive entropy production, i.e. violation 
of entropy conditions for ©. We are currently able to treat the functional 
of [H |37] but not the weaker one of [281 HH] • For the latter we would need a 
single-entropy version of the structure results of [18J HO] for convex conser- 
vation laws. Therefore our result in its present form is not sufficient to deal 
with ASEP large deviations, but presumably effective for © in the limit of 
[37] . We hope it could be used also for vanishing viscosity WASEP ([27]), 
which is formally similar to ©. The boundary terms are a generalization 
of those introduced in [13] for the ASEP case (fT5"]) . and measure violation of 
BLN boundary conditions ©. Although we only study the quasi-potential 
for fluces satisfying ©, the dynamic action potential that we define is rel- 
evant to general fluces. In particular, like the bulk functional of [4J, our 
definition of the boundary part does not assume any convexity for /. This 
general joint bulk/boundary action functional has a natural variational for- 
mula. The quasi-potential on the torus for the inside part of the functional 
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was derived in [5] and yields the equilibrium entropy (I13p , again in agreement 
with the uncorrelated steady states of ASEP-like systems on the torus. 

In Theorem 12.21 we obtain a closed variational expression for the nonequilib- 
rium quasi-potential relative to ©-(EJ. To do so we need a joint symmetry 
assumption on / and h, see (j3"3"]) -(]4"2 ]) . satisfied in particular by ffTo"j) . This 
assumption can be interpreted as the flux being an even function of the 
chemical potential 9 = h'(p) (that is the variable dual to p w.r.t. the en- 
tropy h), see ( H3|) . In a similar spirit, it was observed in [10] that some 
(more stringent) relation between diffusion and mobility was required to get 
a closed expression for the quasi-potential in (15|)- (ITU1) - (ITT|) . Our expression 
is a natural generalization of |23j, which gives a candidate static LD func- 
tional for ASEP-like systems satisfying the symmetry assumption. Besides, 
by definition of the quasi-potential, for given flux /, it yields a whole family 
of Lyapunov functionals for the evolution problem ([2]) (|3]) (namely for all 
entropies h compatible with (j42p ). The adjoint evolution problem is then 
described in Theorems I2.3H2.4I for the space-time reversal of the fluctuation 
path p(., .). Space-time reversal is more natural here than time reversal alone, 
because it leaves fl2]) invariant, hence relaxation paths can be described with 
reference to the original problem (j2]). Our results also cover the phase tran- 
sition line pi < p r , f(pi) = f(p r ), where there is a continuum of stationary 
solutions for ©-Q parametrized by y G [0, 1]. 

There is a connection between our results and those of [8] [9]. For given 
viscosity v in fl9l- f|T0l) - f|Tll . in the case (fl2l) . a quasi-potential V v is obtained 
in [8] [9] from the diffusion- like action functional I v . For (j2])-([3]), we obtain 
a quasi-potential V from the singular action functional /. Despite the dif- 
ference in nature between I u and /, there is strong evidence that I v — > I 
holds in the sense of T- convergence. T-convergence of the inside part of the 
functional is established on a set of "nice" functions (|4J). On the other 
hand, T-convergence of the boundary part is suggested in [13] . With a result 
similar to [I] at the boundary, for (|T5|) and the action functional I restricted 
to nice functions, one could deduce Theorem 12.21 from [5] [9] in the case 
ffT2^ - ([T5]) . Conversely, Theorems I2.3f[2~4l would determine the limit v — > 
of the fluctuation paths in [8] [9]. Strictly speaking this is true outside the 
phase transition, which is not seen at positive viscosity (but can be recovered 
differently by sending boundary conditions to ±oo, see [TT]). 
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We end up with an informal description of the fluctuation paths described 
in Theorems I2.3H2.41 A significant difference with known results at positive 
viscosity, and with the equilibrium case of [5J is that, outside the maximum 
current phase, the fluctuation is reached in finite time. In the regime p\ > p r , 
the quasi-potential is expressed as a supremum over an auxiliary function 
that depends nonlocally on p(.). First, the reversed evolution of this auxiliary 
function is determined as the entropy solution of §2§ with space-reversed BLN 
conditions ()3]). Next, the reversed (unique) fluctuation path is obtained as 
the entropy solution to <^ with space-reversed, time-dependent BLN bound- 
ary data depending on the boundary values of the auxiliary function. The 
first step is reminiscent of [S] [B], but the specific second step is necessary be- 
cause (unlike at positive viscosity) one cannot invert the relation p(.) i— > F(.). 
In the regime pi < p r , the quasi-potential is expressed as an infimum over 
an auxiliary position y 6 [0, 1]. We show that to each minimizing y one can 
associate a unique fluctuation path, and that there are no other fluctuations 
paths than these. This non-uniqueness phenomenon can be understood as 
the v — Y limit of the one exhibited in [9]. The adjoint evolution problem for 
the reversed path here is quite different from prior results and consists of two 
distinct phases. In a first step an antishock <p(pi) > f(p r ) is instantly cre- 
ated at y and travels up to one of the boundaries at Rankine-Hugoniot speed. 
During this step the solution is entropic outside this antishock and follows 
space- reversed BLN boundary conditions pi,p r at each boundary. Once the 
antishock reaches a boundary there is an abrupt change in the rule. Suppose 
e.g. f(pi) < f(p r ), so the stationary state is pi. Then the solution becomes 
entropic inside and the BLN data p r at one boundary changes to pi, so that 
pi is now imposed on both sides and the solution eventually relaxes to pi. On 
the phase transition line f(pi) = f{p r ) the antishick does not move, and the 
stationary state eventually reached is the shock (pi,p r ) at initial minimizing 
position y. 

The paper is organized as follows. In Section [21 we define the setting and 
state the main results. Proofs are carried out in Sections [3] and IU respec- 
tively for the shock regime pi < p r and the rarefaction regime pi > p r . These 
two sections are similarly structured into three subsections, that respectively 
contain preliminary computations and each of the two inequalities between 
the quasi-potential and the static functional. In appendix [AJ we establish 
a variational formula for the dynamic functional, which implies it is a good 
functional in the terminology of LD theory. Eventually, in Appendices [B] 
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and O we explicitely compute fluctuation paths leading to arbitrary uni- 
form profiles, which involves interaction of Riemann waves issued from the 
boundaries. 

2 The setting and results 

Let (a, b) be a time interval, — oo < a < b < +00. We define an action 
functional I( a ,b) on the set L°° ,K ((a,b) x (0,1)) of [0, K]- valued Borel func- 
tions p(., .) on (a, b) x (0, 1). This functional consists of a bulk term 1° and 
boundary terms I 1 and I r : 

I(a,b) = l\ a ,b) + I{a,b) + ^(0,6) 

For its definition and basic properties, we assume first that / in ([2]) is a func- 
tion of class C 1 that is not affine on any nonempty open subinterval of [0, K]. 

Dynamical functional: bulk part. An entropy is a convex function 
T] : [0, K] — > R, and the associated entropy flux q is defined up to an ad- 
ditive constant by 

q'(p) = V'(p)f(p) (17) 

(r),q) is called an entropy-flux pair. The Kruzkov entropy-flux pair ([22]) is 
given by 

Vv (p):=(p-v) + , q v (p) := l p>v [f(p) ~ /(«)] (18) 

Let fl be an open subset of M 2 . The //-entropy production of p(., .) G L oc,,is: (r2) 
is the distribution p, v [p(-, .)] on ^ defined by 

/i,[p(.,.)] :=d t ri(p(t,x)) + d x q(p{t,x)) (19) 

In particular, (id, f) is an entropy-flux pair, and /?(.,.) is a weak solution 
to (J2J) in f2 iff. j^id [/?(., •)] = 0. Notice the behavior of p v under space-time 
reversal: 

Pr, [p o§\ = -pr,[p\ o $, $ : (t, x) ^ (-t, 1 - ar) (20) 

Let M(Q) denote the set of Radon measures p, on fl that are locally bounded 
up to the boundary of Q, i.e. fl f2) < +oo for every compact subset 

K of M 2 . We extend the definition of jl] to the boundary, by calling p(., .) G 
L°° ,K (Q) an entropy- measure solution of 02]), iff. it is a weak solution of ([2]) 
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in Q such that p v [p(-, .)] £ M(fi) for every C 2 entropy rj. The set of entropy- 
measure solutions is denoted by S(Cl). By a straightforward extension of 
(|H Proposition 2.3]), p(., .) G £(fi) iff. there exists a bounded measurable 
mapping v !->■ m p r^(v;dt,dx) from [0, if] to M(fl) such that, for every C 2 
entropy 17, 

,* 

/i^[p(., .)](dt, dx) = / r)"(v)m p (^(v;dt,dx)dv (21) 
io 

m p (. i .) is the kinetic entropy defect measure of p(., .) in the sense of [31] ■ In 
the sequel, (h, g) denotes a given entropy-flux pair with uniformly convex 
h G C 2 ([0,K]), i.e. h" > c > 0. For notational simplicity, we often write p v 
for /x,j[p(., •)], p for /i/j, and m for rn p (\. Now assume f2 = (a, 6) x (0, 1). We 
set 

r +00 if P (.,.) ^£((0,6) x (0,1)) 

I(a,b)(p(; ■)) = \ I I h"(v)m + (v;dt,dx)dv if p(., .) G £{{a, b) X (0, 1)) 

L J [0,isT] J(a,b)x(0,l) 

(22) 

where m + (t> ; (it, efcc) denotes the positive part of m(v; dt, dx). Thus, zeroes of 
l9 a b -K are exactly entropy solution to (JSJ) in (a, 6) x (0, 1). This functional was 
introduced in [3] and shown in |37J to be the LD rate function for a stochastic 
perturbation of @. 

Dynamical functional: boundary part. We set 

Il >b) [p(.,.)} = [ i l (p(t,0 + ),Pl)dt (23) 

J a 

W(">-)] = f i r (p(t,n,p r )dt (24) 

a 

with functions ^ and i r defined below. We point out that f l23l) - fT24l) makes 
sense for p(.,.) G £((a,b) x (0,1)). Indeed, £(fi) is contained in the set 
Q(fi) of quasi-solutions of (T5]) introduced in [39] EQ], i.e. functions p(., .) G 
L°°^ K {yi) such that /i^[p(., .)] G M(ft) for a dense subset of v G [0,K\. For 
such functions, existence of boundary traces in L 1 sense (see (186]) below) is 
established in [3D]. Given an entropy- flux pair (17, g) and a continuity point 
Po G [0, K] of 77, we define the relative entropy-flux pair by 

v(P,Po) '■= v(j>) -v(Po) - v'(Po)(p-Po) (25) 
g(p,Po) := q(p)-q(po)-v'(Po)[f(p)-f(po)] (26) 
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In particular, for (r] v ,q v ) in (IT5|) . T] v (p, p ) and q v (p, Po) are we U defined for 
f 7^ po, and, for every C 2 entropy-flux pair (77, g), 

/■X pK 

v(p,Po) = v"( v )Vv(p,Po)dv, q(p,p )= r)"(v)q v (p, p )dv (27) 
Jo Jo 

We then set in (|23l)-(|2lD 

i{p,Pi) ■■= / h"(v)max[q v {p,pi),0}dv (28) 
./o 

i r {p,p r ) ■= I h (v)max[-q v (p,p r ),0}dv (29) 

Define the sets 







£ + { Pl ) = {p€[0,K\: i l (p,pi) = 0} (30) 
£-{p r ) = {pe[0,K\: f(p,p r ) = 0} (31) 

It follows from (1271) that p G £ + (pi) (resp. p G £~(p r )) iff. q(p, pi) < (resp. 
—q{p,p r ) < 0) for every entropy-flux pair (r/,g). Thus (see e.g. (HI Chapter 
15]) p(t,0 + ) G £ + (pi), p(t,l~) G £~(p r ) are equivalent to BLN boundary 
conditions for ([3]). 

In Appendix El we establish a variational expression for I( a ,b) and prove that 
I(a, t b) is l.s.c. with compact level sets with respect to the local Z^-topology 
on L°° ,K ((a, b) x (0, 1)). One can show (see [10] and beginning of Subsection 
I3TTD that {If ab) < +00} is contained in C°([a, b] n R, L^O, 1)) . 

Stationary set. The set S of stationary solutions of is studied 

in [2] and the following results established. Let M. denote the set of min- 
imizers (resp. maximizers) of / on [pi,p r ] (resp. [p r ,pi]) if p\ < p r (resp. 
Pi > Pr)- Then we have 

Theorem 2.1 ([2]) (i) If pi < p r (resp. pi > p r ), S is the set of piecewise 
constant, nondecreasing (resp. nonin creasing) , Jv[ -valued functions p s {.) on 
[0,1]. In particular, if M. is a singleton, S is a singleton whose unique 
element is the function p s (.) with constant value 

R(p h p r ):=h Igmhl ^ f f % Pl i Pr (32) 
[ argmax [ft pA f if p t > p r 
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(ii) If M. is a singleton, the entropy solution p(t, .) to (d))-([3|) with Cauchy 
datum po(-) L°° ,K (0, 1) converges to p s (.) in L 1 (0, 1) as t — > oo. // in 
addition R(pi, p r ) is not a local extremum of f , there exists T e [0, +oo) 
such that p(t, .) = p s (.) /or a// 1 >T. 

Exclusion-like case. From now on we assume that / is an "exclusion-like" 
flux function, i.e. satisfies ([1]). Thus there exists a unique decreasing function 
<p : [0, K\ ->■ [0, K\ such that 

f°<P = f (33) 
In this case, by elementary computations, one finds the following expressions 
for flU-flU}, (PP-flEP and S. 



if p = pi or <p(p*) < p < K 

g{p, Pl ) if 0<p<p* (34) 

9(<P(p),Pi) if P*<P<^(Pi) 



if p* < p < K 

g(p, Pl ) if 0<p<<p(p z ) (35) 

g(p) - if y(pj) <p<p* 



if < p < p* 

-g(p,p r ) if f(p r )<P<K (36) 

«(^(p)) - O(p) if P* < P< ip(Pr) 



if p = p r or < p < <p(p r ) 

-o(p,p r ) if p*<p<l (37) 

-9{<P(ft),Pr) if V(Pr) < P< P* 

|)- (|37]l were introduced in [T3] in the special case (IT5"|) . For sets of admis- 
sible boundary values, we obtain 

vw \ [p*,K] if Pl >p* v 7 
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E (p r 



(39) 



[0,p*\ if Pr <p* 

{ Pr }U[0,p*] if Pr >p* 

Finally, the stationary set S reduces to the four phases introduced in [21]: 



S 



' {p s (.) = pi} 

{Ps(-) = Pr} 
{Ps(-)=P*} 



if pi < p r and f(pi) < f(p r ) 

or p* > pi> p r and f(pi) > f(p r ) 

if pi < p r and f(pi) > f(p r ) 

or pi> p r > p* and f(p r ) > f(pi) 

if Pi> P* > Pr 



{ {#(.), y e [0, 1]} if p x < P r and f( Pl ) = f(p r 



(40) 



where 

#(.)=Pll(0,y)+Prl(y,l), V & [0, 1] (41) 

The four cases in (14"0"|) are respectively ( [2 lj ) the low-density (LD), high- 
density (HD) and maximal current (MC) phase, and the phase transition 
line between LD and HD. 



Static functional. We now make the symmetry assumption (see (1531 ) 



h'\p) = -h"{ V {p)) V '{p) (42) 

Since (142]) and the definition of I{ a ,b) is unchanged by adding a linear term 
to the entropy h, we may assume without loss of generality that h = h(p, p*) 
(see ([25])) is the relative entropy with respect to p*. With this choice (1421 
is equivalent to h'((p(p)) = —h'(p). Let 9 = h'(p) be the chemical potential. 
The Legendre-Fenchel dual L(8) = h*(8) can be interpreted microcsopically 
as the moment generating function under the equilibrium invariant measure 
of an underlying particle system with equilibrium entropy h. The latter can 
be thought of as the equilibrium LD rate function for the empirical particle 
density. Let j{9) = f{h'^ 1 {9)) be the flux as a function of the chemical 
potential. Then (1421 is equivalent to 



K-o)=m 



(43) 



Let 



k(p) := / (p(r)h"(r)dr, K(p) := h(p) - ph'(p) + k(p) (44) 
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It follows from (@2J) that 



k(p) = <p{p)h!{p) + %(p)], K{p) = [<p(p) - p\h'{p) + h(p) + %(p)] (45) 
and 

K(<p(p)) = K{p) (46) 

K is increasing on [0, p*} and decreasing on [p* , K], and there exists a unique 
function L : [0, f(p*)} — >■ K. such that X(p) = L(f(p)). L is an increasing 
function. We now define on L°°'^((0, 1)) a generalization £[.] of the station- 
ary LD functional introduced in 



First case: pi < p r . For p{.) G L°°> K ((0, 1)), and y G [0, 1], set 

%(.),!/] := / [%(x))-p(x)/i'(pO+^(p/)]^+ / [%(x))-p(x)fc'(p r )+fc(Pr)]<fo; 

(47) 

%(.)] := inf %(.),!/] -min(ir( A ),^(Pr)) (48) 
!/e[o,i] 

We denote by 3^[p(-)] the se ^ °f minimizers y in (148]) . that is the set of 
minimizers y of y i— >■ p(z)dz — p c y, where p c is given by 

_ . fc(p r )-fe(pz) (49) 



rt ' /i'(p r ) - h'(pi) 
Second case: p x > p r . For F(.) G L°°' x ((0, 1)), set 

%(.),F(.)] := [ [h(p(x))-p(x)h'{F(x))+k(F{x))]dx (50) 
Jo 

%]:= sup S[p{.),F(.)]- sup K(p) (51) 

where J-"(pz,p r ) is the set of nonincreasing functions F on [0,1] such that 
Pi > F(Q+) > p r > F(l~). As in [23j, the unique minimizer in (1511) . which 
we will denote by F p ^, can be described explicitely by an envelope construc- 
tion, see beginning of Subsection 14.31 One shows as in [23] that 5[p(.)] = 
for p(.) G S, and 5*[p(.)] > for p(.) G" S. It is not difficult to see that 
the functionals (j4"8]) - (l5"Tj) are continuous with respect to L l ((0, 1)) norm and 
l.s.c. in weak* toplogy a(L°°' K (0, 1), L 1 ((0, 1)). 
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Quasi-potential. Let lZ\p{.)} denote the set of functions p(., .) G L°°' K ((— oo, 0) x 
(0,1)) such that lim^o p(t, .) = p(.), and lim^_ 00 p(t, .) = p s (.) G S in 
L 1 ((0, 1)). For T G (— oo, 0), let 1Zt[p(-)] denote the set of functions p(., .) G 
L°°^((-oo,0) x (0, 1)) such that lim t ^ p0, •) = p(-) in ^((0, 1)), and there 
exists p s (.) G S for which p(T, .) = p s (.) for all t < T. Define 

p(.,.)eft [/>(.)] 

^ T [p(.)] := inf rn /(-oo,o)b(.,)] (53) 

We can now state our main results. The first theorem identifies the quasi- 
potential with the static functional, and states that it is achieved in finite 
time outside the MC phase. The next two theorems characterize optimal 
paths for the quasi-potential. 

Theorem 2.2 For every p(.) G L°°' x ((0, 1)), 

%(•)] = V[p(.)] (54) 

Outside the MC phase, we also have 

S[p(.)} = mmV T [p(.)} (55) 

We denote by v(p~,p + ) the Rankine-Hugoniot or velocity of a discontinuity 
(p~ 7^ p + ) in ([2]), given by the flux continuity relation 

f(p + )-vp + = f(p-)-vp- (56) 



Theorem 2.3 Assume p\ < p r . Then there exists a mapping 

(p(.) G L°°' K ((0, 1)), y G y[p~(.)]) H- p(., .) = M*\p(.)] G L^((0, +oo) x (0, 1)) 

with the following properties, where y[p{-)} denotes the set of maximizers of 
V ^ Jo P( z ) dz - PcV- 

(0) M^[p(.)] ^ iVF 2 [p(.)] for two distinct y\y 2 G y[p(.)). 
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(1) p(.,.) G ^'^((-00,0) x (0,1)) achieves equality in and (5^) iff. 
it is of the form 

p(t,x)=p(-t,l-x) (57) 

where p(., .) = M x ~ v [p(.)] for some y G 3^[p(-)L P(-) = P(l — •)■ 27jen 
we nave 3^[p(t, .) = {yt} for all t < 0, where y% — 1 — y-t, yo = 2/ ; with y, 
given by 155]) below. 



(2) M*[p(., .)] zs the unique weak solution of m £((0, +oo) x (0, 1)) that 
satisfies the following conditions (a)-(c). Besides, y[p(t, .)] = {y t } holds for 
all t > 0, with y t given by / fjgj) . 

(a) p(t, .) -> p(.) m ^((0, 1)) as t -> 0. 

(b) Let r = {(t,y t ) : t e (0, ^)} ; wnere 

y t :=min(l,(y + <)+), t G [0, +oo) (58) 

v := v(<p{pi),(p{p r )) (59) 

0* : = max(—y/v, (1 — y)/v) G [0, +oo] (60) 

rai/i 0/0 = 0, is the time fory, to reach a boundary. Then p(., .) is an entropy 
solution to (HP in [(0,9$) x (0, l)]\r, and satisfies 

p(t,0 + )eS + ( Pr ), p(t,l-)eE-(jH) (61) 

p(t, y t -0) = <p{pi) > p(t, y t + 0) = <p{p r ) (62) 

(c) If 6$ < +oo, i.e. f(pi) 7^ f(Pr) oryE {0, 1}, p(., .) zs an entropy solution 
to (U|) m (6^, +oo) x (0,1), and satisfies boundary conditions 

p(t, y a J = Pb, p(t, (1 - yoo)- a ) g £-°( Pb ) (63) 

~ ~ _ / if f(pt) < f(p r ) or f( Pl ) = f(p r ), y = 

2/00 Ve " I 1 /(A) > /(Pr) ^ f( Pl ) = f(p r ), y = 1 [M) 

is the boundary reached by y, at time 6^ , a = + and —a = — if y^ = 0, 
a = — and —a = + if y^ = 1, p& = p^ m £/ie /irst case of (64\ ), and p b = p r 
in the second case of 
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Theorem 2.4 Assume pi > p r . Then there is a unique p(., .) G L OD,K ((—oo,0)x 
(0, 1)) that achieves equality in [blty and ( [5^] ), which is defined as follows. 
Let 

G(Q,x) = V [F p{ . ) (l-x)} (65) 
Define G(., .) on (0, +oo) xl as the entropy solution to (TJ)) with initial datum 

(5(0, .) = V2(p r )l(-oo,o) + G(0, .) l(o,i) + y(^)l(i,+oo) (66) 

T/ie restriction ofG(.,.) to (0, +oo) x (0,1) is the entropy solution to (TJ|) 
with initial datum G(0, .) and BLN boundary conditions 

G(t, 0+) G S + (if( Pr )), G(t, 1-) G (67) 



T/ien p(.,.) is given 61/ ( [571 ), where p(t,x) is the unique entropy solution to 
([IP with initial datum 

p(O,.) = p(l-0 (68) 

and BLN boundary data 

pit, 0+) G £+[<p(G(t, 0+))], p(t, 1-) G £-[^(G(t, 1-))] (69) 

Besides, ([22P remains true at later times, i.e.: 

G(t, x) = ^[F p( _ til _ (l - a;)], Vt > 0, a; G (0, 1) (70) 

A functional closely related to (}2"2"]) was introduced in [281 SS] as the LD rate 
function for TASEP on Z. Namely, 

70 r n( \] = l+°° if &£h((a,b) x (0, 1)) 

^)IPU-JJ I M + ((a)6) x (0jl)) jf pG;0 G5h((aj6) x (o,i)) ^ 



J ( %>J ( % on ^((a, 6) x (0,1)) (72) 



where £/i(f2) is the set of weak solutions of (j2J) such that p/i lies in M(fi). 
By convexity, 

l (a,b) — J (a,b) 

If / is strictly convex or concave, the result of (18] implies that zeroes of J? a ^ 
are exactly entropy solution to 02]) in (a, b) x (0, 1). For such / one has 

J { % = /(% on BV loc ((a, b) x (0, 1)) (73) 

where BVi oc (fl) denotes the set of functions p(., .) G L°°(f2) such that d t p 
and c^p lie in M(Q), the set of Radon measures on Q. If / is neither 
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strictly convex nor strictly concave, 1° and J° are quite different, as a 
single entropy is not enough to characterize entropy solutions ([31]). In 
fact we will prove our results more generally for a bulk functional 19 b ^ : 
L°°' K ((a, 6) x (0,1)) — > [0, +oo] which enjoys the following properties: (i) 



J(a, b ) < ^ L°°' K ((a,b) x (0, 1)); (ii) I? ab) = J^ b) on BV loc ((a, b) x (0, 1)); 
(hi) (J°)^ ) ([0,+oo))cf((a,6)x(0,l)). 



We will denote by J( a ,b) the functional obtained by replacing 1° with J° 
in (fT6j) . Since we do not know if J° satisfies assumption iii), in view of 
flH])-(|2l]), J {a<b) is not defined a priori on L°°' K ((a,b) x (0, 1)), but only on 
£((a,b) x (0, 1)). Our results would be valid for J( a ,b) if condition (iii) could 
be proved in this case. The result of [T5| suggests the possible conjecture 
Sh((a, b) x (0, 1)) = £((a, b) x (0, 1)) for any strictly convex entropy h, which 
would imply iii) for J( a ,b)- 

3 The case p\ < p r 

3.1 Analysis of entropy production 

Let p(.,.) G S((a,b) x (0,1)). It will be convenient to extend it by giving 
it constant uniform values po,Pi, respectively for x < and x > 1. This 
extension is in general no longer a weak solution of ([2]). However, by [T5| 
Proposition 2.1], it still lies in the set Q((a,b) x R) of generalized entropy 
solutions in the sense of pEH], i.e. functions p(.,.) e L°°' K ((a,b) x R) such 
that p v [p{-i •)] £ M((a, b) x R) for every C 2 entropy 77. 

Recall that a bounded Borel function F on M. d is said to have an approx- 
imate limit l n F(x) at x G M. d along the normal n G R d , ||n|| = 1, if 



where B(x,e) denotes the open ball of radius e centered at x and H n (x) = 
{y G R d : n.(y — x) > 0} is the normal hyperplane containing x and n. x is 
called a jump point of F if, for some n, approximate limits l± n F(x) exist and 
are distinct. One can show that the pair {n, —n} is uniquely determined, x 




(74) 
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is called a point of vanishing mean oscillation (VMO) of F if 

lim e~ 2 [ \F(y) - F(x, e) \dy = (75) 

JB(x,e) 

where F(x,e) is the mean value of F on B(x,e). 

Let (77, g) be a C 2 entropy-flux pair. By [19], /?(•,.) 6 (?((a, 6) x M) has 
the following property: there are subsets J and V of (a, 6) x R such that the 
complement of J U V has ^-measure 0, J is rectifiable with a local normal 
n = (n t ,n x ), p(., .) has approximate limits p + 7^ p~ on J (p + being the limit 
in the positive n-half-space), and p(., .) has VMO on V. The trace of entropy 
production on the jump set J is 

p r ,[p(-, -)}\-J = 7r v (p~,p + , n)U l ^J (76) 

where 

7r v (p~,p + ,n) = n x [q(p + ) - q(p~)\ - n t [r](p + ) -T](p~)\ 

For (r],q) = {id J), © yields n x + 0, % x -a.e. on J n {(-00, 0) x (0,1)}, 
and the Rankine-Hugoniot relation (I56p for the local speed v := —n t /n x . By 
convention we orient the normal so that n x > 0. Note that n x 7^ implies 
no jump along time lines. Together with [39] this implies £((a, b) x (0, 1)) C 
C°([a, b], -t/ 1 ((0, 1)). We now focus on (77, q) = (h, g) and set 

n(p-,p + ) = [g(p + )-v(p-,p + )p + }-[g(p-)-v(p-,p + )p-] (77) 

h"(r)[f(r) - f(p+) - v(p-,p + )(r - p+)) (78) 



Let 



a(p,0:= h"(r)[f(tp(r))-f(p)-av(r)-p)]dr (79) 
J pi 

Notice that u(p~ , £) = o~(p + , £) if C — ^(p - ! P + )- 

Lemma 3.1 Let p(., .) G £((a,b) x (0,1)), and T = : £ G (a, b)}, 

where y. : [a, 6] — > [0, 1] zs a Lipschitz trajectory. Set S t = S[p(t, -),yt\- Then 

s b -s a - J (0 ,6)[p(., .)] = / - p~ [M) x (0, i)\r] 

l(o,i) (Vt)n(p(t, y t -),p(t, y t +))~dt (80) 



18 



where 



d(t) 
d\t) 



d°(t) + dd(t) 
^(o,i)(yt)cr(p(t,y t ±),y t ) 



(81) 
(82) 



and dd(t) is given respectively when y t G (0, 1), y t = 0, yt = 1, by 



[g(p(t, + ); Pl ) - i\p(t, 0+); P i)\ + [-g(p(t, 1"); p r ) - f (p(t, 1"); p r )] (83) 

[<?(p(t, + ); p r ) - i l (p(t, 0+); p,)] + h^(p(*, I")! - i r (p(t, I"); Pr)] (84) 
[S(p(t, + ); pO - i l (p(t, 0+); pO] + [-g(p(t, 1"); p,) - f (p(i, 1"); p r )] (85) 



Remark. The above limits can be understood as approximate limits provided 
the last integral in ( 180|) is restricted to T D J. Due to the extension of p(., .), 
VMO points at the boundaries have approximate limit inside, and jump 
points at the boundaries have a normal parallel to the space axis. Hence 
approximate spatial limits are defined a.e. at the boundaries. Alternatively, 
since p(., .) G Q((a, b) x (0, 1)), the limits in Lemma I3TT1 exist in the L 1 sense 



my- 

esslim / \p(t, e)-p(t, + )\dt = ess lim / \p(t,y t +e)-p( y t,y t +0)\dt = (86) 



for every bounded interval / C (a, b). 

We need a version of Lemma 13.11 for a mollified version of S[.}. Let ip : 
[0, 1] — > [0, 1] be a nonincreasing function such that -?/>(0) = 1 and ^(l) = 0. 
We set 



We denote by BVi oc (Q) the set of functions p(.,.) G L°°(Q) such that d t p 
and d x p lie in M(Q). 
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Lemma 3.2 Let ip G BVi oc ((a, b) x (0, 1)) . Assume ifj(t, 0) = 1 andip(t } 1) = 
for a.e. t G (a, b), d x ip < 0, \d t i[)\ <C —d x i/;, and d x ip(dt,dx) = d x ip(t,dx)dt. 
Let S t = S^ t ,.)[p{t, •)]■ Then 

S b -S a -J {a>b) [p{.,.)}= [ d4,{t)dt-fj,-[(a,b) x (0,1)] (89) 

J a 

where (with dd given by §EB\) ) 

d^{t) := d^it) + dd(t) (90) 

d;(t) := J\(p{t,x) t ~£$l'*f ) i-dMt,dx) (91) 

The next lemma will be combined with f[8"0~j) - fl8"9"j) to show that dd(t), d°(t) 
and d^(t) are nonpositive for suitable choices of y, and if), and derive condi- 
tions from these quantities being equal to 0. 

Lemma 3.3 Let p~ , p + G [0, if]. Then: 

(0) p~ < p + (resp. >) iff. n(p~,p + ) < (resp. >). 

(1) (a) Assuming p~ < p + , 

a(p ± ,v(p-,p + )) + 7T(p-,p + )<0 (92) 

with equality iff. p~ = (p(p r ) and p + = <f(pi). (b) 

sup a(pj'(p))<0 (93) 
pe[o,K] 

(ii) i\p,pi) > g(p,pi), with equality iff. p G £~(pi); i r (p,p r ) > -g(p,p r ), 
with equality iff. p G £ + {p r ). 

(Hi) Assume p < p* . Then —g(p,pi) < i r (p,p r ), and 

- 9(P, Pi) = i r (p, Pr) f(pi) < f(pr) and p = pi (94) 
Under P^] ) ; we have —g(p,pi) = i r (p,p r ) = . 

(iv) Assume p > p* . Then g(p,p r ) < i l (p,pi), and 

g(p, p r ) = i\p, pi) O f(pi) > f(pr) and p = p r (95) 
Under ( GSP , we have g(p,p r ) = i l (p,pi) = 0. 
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In the proof of Lemma 13.31 and various other places we use the following 
lemma, which is immediate from assumptions ([T]) on /, (j!7p and (1331) . 



Lemma 3.4 (i) Let p < p* . Then g(p,p ) > for p < p* , p ^ p ; and 
g(p, p ) < for p> ip(po), p p . 

(li) Let p > p*. Then g(p,p ) < for p > p*, p ^ p ; and g{p,p ) > 
for p< <p(po), p ^ po- 

Proof of lemma [X71 Apply the generalized Green's formula Q15J) to the 
divergence-measure fields (h(p), g(p)) on (a, b) x (0, 1), and (p, /(p)) on either 
side of T. For the latter, use (j2J) and note that y t = a.e. on y _1 (0) Uy~ 1 (l). 
The result then follows from simple computations using (IT7|) . fl33|) . (jUj) and 
fl56|) . We eventually use f!76|) to decompose p~((a,b) x (0,1)) arising from 
the first Green's formula into the second and third term on the r.h.s. of (180]) . 



Proof of lemma \3.2i The proof is similar to that of Lemma I3.1[ but here 
the second Green's formula must be applied to if)(t, x)[p(t, x), f(p(t, x))] on 



Proof of lemma \3.3\ 

Proof of (o). The claim follows from flTHj) and strict convexity (resp. concav- 
ity) of h (resp. /). 

Proof of (i). (b) Follows from strict concavity of / and strict convexity 
of h. For (a), setting p = <p(r) in (179!) and and using (|42j) . we obtain, with 



p i— > /(p) — £p is a strictly concave function that vanishes for p = p ± , because 
£ = v(p~,p + ). Hence the integrand is negative for p ^ [p~,p + ], and positive 
for p G (p~, p + ). This yields f[9"2"j) . and shows that equality occurs if and only 



□ 



(a,b) x (0,1). 



□ 



i = v(p ,p + ), 
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if Lp(p r ) = p and <p(p{) = p + . 

Proof of (ii). We establish the result for i l , the proof for i r being simi- 
lar. We distinguish two cases: 

First case: pi < p* . Then by (pUD and (|3"0]) . we have g(p,pi) = i l (p,pi) if 
p < p*, that is p G £~(pi). When p > <p(p;) and p ^ pi, f[3~4"l) and Lemma [BT^I 
imply g(p,pi) < = i l (p,pi). For p* < p < <p(pz) we have, by ([26]) and (|33|) . 
i l (p,pi) —g(p,pi) = g((p(p)) — g(p) > 0, where the inequality follows from (o) 
with (r],q) = (h,g). 

Second case: pi > p*. By (]3"5T) and ([3"U]) . we have g(p, pi) = i l (p, pi) if p = pi or 
P < viPi), that is p G £~(pi). By fl35|) and Lemma I3T4"] g(p,pi) < = i l (p,pi) 
if p > p* and p ^ p*. For <p(p,) < p < p*, ([34]), © (03]) and Lemma [331 
imply z'(p, p^ - g(p, p t ) = g(<p(p),pi) > 0. 

Proof of (in). We distinguish two cases: 

First case: p\ < p*. For p < p*, by Lemma [3 A\ 

-g(p,pi)<0<i r (p,p r ) (96) 

This is an equality for p = pi and f(pi) < f{p r )- Indeed, the latter means 
either pi < p r < p* , or p r > p* and pi < (p(p r ). Then, by (!36|) - (l37|) and 
Lemma I3.4[ 

-g(Pl,Pl) = = i r (pi,Pr) 
Conversely, assume —g(p,pi) = i r (p,p r ). By f[9"6"]) . we must have 

g(j>,Pi) = = i r (p,p r ) (97) 

Since p < p*, by Lemma E3, the first equality in f[9~T[) implies p — pi- Then, 
again by Lemma I3T4] the second equality i r (pi,p r ) = in (|97|) holds iff 

f(Pl) < f(Pr). 

Second case: pi > p*. If p < <p(pz), then Lemma implies —g(p,pi) < 
< i r (p,p r ). Assume now <p(p«) < p < p*. Then 

#(p, Pr) - ff(p, P«) = h'(p r )f(p r ) - g(p r ) - [h'(pi)f(pi) - f{pi)\ 
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(h'( Pr ) - h'( Pl ))f( P ) 

= r h"(r)(f(r)-f(p))dr 
J pi 

Since <p(pz) < p < p* and r G [pz,Pr], the integrand is negative. Hence 

-g(p>pi) < -g(p,Pr) < i r (p,Pr)- 

The proof of (iv) is symmetric to that of (iii) and thus omitted. □ 

3.2 S < V and characterization of minimizers 

In this subsection, we establish S < V in (I53|) . and prove that every p(., .) 
that achieves equality in fl52|) and fl5^|) necessarily satisfies conditions of The- 
orem 12.31 In the next subsection we shall prove actual existence of such 
p(.,.). Before going into details, we explain the main idea. Let p(., .) G 
£((-oo,0), x(0, 1)) with lim t ^_oo p(t, .) = p s (.) G 5, p(0, .) = p(.), and 
{?/*} = y[Ps(-)}, that is, y* = 1 if p s (.) = p h y* = if p s (.) = p r , = y in 
the case (Hill . Let (y t , t < 0) be a [0, l]-valued a Lipschitz path such that 
lim^oo^ = y*. Then 

%(.)] = S[p(.)} - S[ Ps (.)} < S[p(0, .),y ] - S[p s (.),y*\ 

If we find a path y. such that d°(t) < and dd (t) < in ( 180|) (with a — >■ — oo 
and b = 0), we obtain the upper bound in f|5^|) . It follows from ( |80|) and 
Lemma [3731 that <i (t) < and dd(t) < if the path y, enjoys the following 
properties: 

(i) Whenever p{t,yf) are different, y t is the Rankine-Hugoniot speed fl56|) . 

(ii) Whenever p{t,yf) are equal, y t is the characteristic speed f'(p(t,y t ±)). 

(iii) p(t,0+) > p* a.e. on {j/ t = 0}, p{t, 0+) < p* a.e. on = 1}. 

(iv) y, sees only shocks, i.e. p(t,y t —) < p(t,y t +) a.e. 

Note that (i) is automatically satisfied by any path y , and (iii) means 
that y. stays at a boundary iff. it is driven against it by the character- 
istic speed. Equality in flHUj) implies that the measure term vanishes and 
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d°(t) = dd(t) = 0. The former gives entropicity of p(., .) outside T. The 
latter yields conditions on p(t,y^), p(t,0 + ), p(t,l~). Using Lemma 13751 one 
sees that these conditions are exactly those of Theorem 12.31 

A path satisfying (i)-(iv) is a (maximal or minimal) generalized forward 
characteristic (g.f.c.) of p(.,.) in a sense slightly wider than [T7] (where no 
boundary and only entropy solutions are considered). Its existence can be 
established with the arguments of [T7J Theorem 3.2] if we know that limits 
p(t,x±) exist in usual sense for a.e. t. Unfortunately, in £((— oo,0] x (0, 1)), 
we only have these limits in the weaker sense (1861) , and we are not able to ob- 
tain existence of the g.f.c. in this case. Therefore the actual proof is a local, 
somewhat more technical version of the above idea, that uses the mollified 
version (1891 of flHDj) . However, for paths that achieve equality in (1521) and 
fl54|) . we can a priori establish existence of genuine limits in certain domains, 
and (l80j) will be used to reach some of the conclusions. 

The first lemma is a simple consequence of the fact that £((— oo, 0) x (0, 1)) 
is contained in C°((— oo, 0], I/ 1 ((0, 1)). Then, all but the existence part of 
Theorem 12.31 is contained in Propositions I3.1H3.3I below. 

Lemma 3.5 Let p(.,.) G £((— oo,0) x (0,1)). Then the following set is 
closed: 

M[p{., .)] := {(t,y) G (-oo,0] x [0, 1] : y G y[p(t, .)]} (98) 
Proposition 3.1 Let — oo < s < t < 0, p(., .) G S ((-oo, 0) x (0, 1)). Then 
S[p(t,.)]-S[p(s,.)} < J (s , t )[p(.,-)] (99) 
In particular, if p(.) G L°°' K ((0, 1)) and p(., .) G TZ[p(.)}, 

%(•)]< J(-oo,o)[p(-,.)] (100) 

Assume (G2|) is an equality. Then for a.e. 9 G (s,t), y[p(9, .)] is reduced 
to a single point y such that (0,y) G J if y G (0, 1), and the following hold, 
respectively if y G (0, 1), y G (0, 1], y G [0, 1), y = l, y = 0: 



p{o,y-) 


= <p(p r ) < p(0,y+) = <p(pi) 


(101) 


p(6,0+) 




(102) 


p(M~) 


G ^ + (p r ) 


(103) 


p(M") 


= Pi, /(PO < /(Pr) 


(104) 


p(M + ) 


= Pr, /(Pi) > /(Pr) 


(105) 
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Besides, 

//-(fi) = (106) 
where Q is the open set defined by 

Q := {(9, y) G (-00, 0) x (0, 1) : y £ y[p(6, .)]} (107) 

Proposition 3.2 Assume p(., .) G IZ[p(.)} achieves equality in |Hj] and (5$ ). 
Then p(.,.) G L 00 '^ ((0, +oo) x (0,1)) given by (57| ) satisfies conditions (a)- 
(c) of Theorem \2. 3\ for some y G 3^[p(.)]> where p(.) = p(l — .). 

Proposition 3.3 Let p G L £X> ' j!r ((0, 1)). Tnen /or every y G !V[p(.)], t/iere 
exists at most one p(., .) G L oo ' f!r ((0, +oo) x (0, 1)) satisfying conditions (a)- 
(c) of Theorem \2.3[ If p(., .) exists, we have y[p(t, .)] = {y t } for every t > 0. 
Besides, p(., .) G L°° ,ir ((-oo, 0) x (0,1)) ^en 6y Hes m TZ T [p(.)} for 
some T < +oo, for which it achieves equality in < f55|) and < f?3)j . If p(.,.) 
exists for two distinct values of y, statement (0) of Theorem \2. 3\ holds for 
these values. 

For the proof of Proposition 13. 2\ we need a local version of [T7J Theorem 3.2] 
to deal with generalized characteristics. We do not repeat the proof, which is 
a straightforward variant of the original one. Let p(., .) G L°° ,K ((— oo, 0) xR). 
Recall that y_ is called a Filippov solution ([25]) of 

yt = f'(p(t,y t )) (108) 

if it satisfies the differential inclusion 

yt G [essliminf /'(p(t,x)),esslimsup/ / (p(t,x))] (109) 

A forward (resp. backward) solution issued from (t ,y ) G (— oo,0) x R 
is a solution w. defined on [to,ti] (resp. [t i , to] ) ? where t < t\ < (resp. 
ti < to), such that ?/ to = y$. Any forward (resp. backward) solution issued 
from (to,yo) lies in a fan between a unique lower forward (resp. backward) 
solution and a unique upper forward (resp. backward) solution. Because /' 
is bounded there is no blowup in fllOSp . hence the upper and lower forward 
(resp. backward) solution are defined on [to,0] (resp. (— oo,t ])- 

Lemma 3.6 Let (to,2/o) G (— oo, 0) x R, and y, be the upper or lower back- 
ward (resp. forward) Filippov solution of U08\) issued from (t ,y ). Suppose 
T is a subset of {— oo,0) such that p(., .) is an entropic (resp. anti-entropic) 
solution of |D) ; i.e. p + [p(.,.)] = (resp. p~[p(.,.)] = 0) in a neighborhood 
°f {(£> lit) '■ t G T}. Then p(t, x + ) = p(t, x~) for a.e. t G T. 
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The underlying heuristics goes as follows. Suppose e.g. p(., .) is anti-entropic 
and y, is the upper forward solution. The only possible discontinuities of p(., .) 
along y, are antishocks, i.e. (since / is strictly concave) p _ = p(t,y^~) > p + = 
pit.yf). Then the Rankine-Hugoniot local speed y t = v(p~,p + ) < f'(p + ). 
Thus by perturbing y. to the right one could obtain a new bigger solution of 
( 11081) that contradicts maximality of y_. 

Remark. Limits in Lemma 13.61 exist in usual sense for all t because ei- 
ther p(t, x) or its space-time reversal is locally an entropy solution and / is 
uniformly concave, thus has locally bounded space variation. 

Proof of proposition VS. 1\ ( llOOp follows from because S[] vanishes on S 
and is continuous w.r.t. L 1 norm. For 9 £ (s,t), let 

D{e) = s[ P (t, .)] - s[ P (e, .)] - j m [ P (; •)] (no) 

and 

A W =1immf P W- J < g + '> (111) 

Let y £ y\p{9,.)\. Then 

Die) - Die + e)< S[p(6 + e, .),y e+£ ] - S[p{9, .),y e ] - J [e ,e + e) [p(, •)] 

for any Lipschitz path |/ on [9, 9 + e] such that yg = y. In particular, taking 
the constant path y=yin f lHUj) . we see that e~ 1 [Di9) — D(Q+e)\ is uniformly 
bounded above. Hence, Proposition 13.11 follows if we show that A(#) < a.e. 
on (s,t), and that A(0) = implies (outside a null subset) y[p(9, .)] = {y}, 
where p(9, .) and y satisfy the conditions of Proposition 13.11 

In the sequel we choose the values po, pi of p(.,.) outside x £ (0,1) such 
that po < P* < Pi- Let J s , resp. J a , denote the shock and antishock set of 
p(., .) in (— oo, 0) x (0, 1), defined as the set of (t, x) £ Jfl ((— oo, 0) x (0, 1)) 
such that pit,x—) < p(t,x+) (resp. >). Let T denote the total subset of 
points 9 £ (s, t) such that 9 is a Lebesgue point of p(6 ) , + ) and p(#, 1~), and 
{9} x (0, 1) C J s U J a U V. In the following we consider (9,y) with 9 £ T 
and y £ !V[p(0, •)]• We set p ± = p(6 l , j/±) whenever (#,?/) £ J s U J a . Let 
£o £ whose value will be chosen below. For 5 > 0, we define functions 
: [0, + e] x [0, 1] -> [0, 1] by ^ ±(5 (^, .) = l [0 , y] and, for u > 9, 

iP ±s (u,x) := * ±s (ix-y)/(u-9)) (112) 
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¥,(0 := i(-oo^)(e) + ^i(eo,eo+<5)(0 



(113) 



:= l(-oo,€ -*)(0 + e ° | - l(g -^o)(0 (H4) 



so that 



■ d x ip± s (u,x) 

d u ip ±s (u,x) 
d x ip ±s (u,x) 



1 :il 



x-y . 
u-9 



L /(£o,£o+<5) 



g - y 
g -y 

u-0 



(115) 
(116) 



with I(a,b) 

1p = V>±5, 



(min(a, 6), max(a, 6)). Since y G !V[p(0, •)] we have, with 



D{9) - D{9 + e)< Sne+e,.)W + e > 01 " ^,.)[p(0, ■)] " <W)[p(- •)] (117) 
Define 



4(«) 



d\(u) 



o (^j5 1(€ ^ 0+5) 



l(€o,«o+<5) 



o (u-e)5 



x-y 
u-9 

x-y 
u-9 



a[p{t,x),Z ]dx (118) 
a[p{t,x),f (p(t,x))]dx (119) 



Since cr(p, £) is uniformly Lipschitz we have, for some constant C = C(f, h) > 
0, 



d°^{u)du — £ 1 d]p(u)du 



< C5 



(120) 



/ d\(u)du < Ca (121) 
Je 



with defined in f l9Tj) . In the following we apply either ( 155]) with ^ = ip±$ 
and £o as specified, or fl80|) with as specified. The constant := l+||/'||oo 
will appear in several places. 

First case. Assume (9, y) G V . Let £ — f'{p)i where p is the mean value of 
p(., .) in the ball Bg jy>e of radius We centered at (9,y). Then 



_-l 



d]p(u)du — e^ 1 / d: 

. +ae Je+ae 



^\u)du 



< CoT x 5- x e- 2 



+£ ry+(£o+6)e 
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\p(u, x) — p\ dxdu 



For small 5 > 0, the domain of the above space-time integral is contained 
in Bq V £ . Thus, by VMO property, this integral vanishes as e — > 0. Letting 
e — > and then (5, a) — >■ (0, 0) in (189]) . using (ii) of Lemma 1331 for the bound- 
ary terms, and (I120p - (I121I) . we find that A(6>) is bounded above by (J9~3"j) . 

Second case. Assume (9,y) G J a U J s . Let £ = v (p,P + ) if ($j 2/) £ ^ 
£o = f'(p + ) if 2/) £ J a - Since <r(p, £) is uniformly Lipschitz, 



Since n 4 = —v(p~ , p + )n x , the domain of the above integral lies in the half- 
plane containing (9, y) and n. This holds also in the case (9, y) G J a , because 
p~ > p + implies £ > v (p~,P + ) by strict concavity of /. Thus, by approx- 
imate limit property, the r.h.s. of f 1 1 2 2 j) vanishes as e — > 0. Thus in the 
case (9,y) G J a we find as previously A(#) bounded above by the l.h.s. of 
( 193|) . Let T' denote the set of 9 G T such that there exists z G y[p(9, .)] for 
which (9, z) G V U J a . In view of the above, to establish A(9) < in remain- 
ing cases, we may assume without loss of generality that 9 G T\T' ■ Once 
the inequality is established for such 9, it will also follow that T' is negligible. 

We now consider (9,y) G J s . We bound the measure term in (1891) as follows. 
Let Z be an arbitrary finite subset of 34* : = y[p(9, .)] fl (0, 1). For each z G Z, 
let n(9,z) denote the local normal at J and v(9,z) = —n t (9,z)/n x (9,z) 
the local velocity. Let H\ and H2 denote nonnegative continuous functions 
on [0, +00) supported and not identically on [0,1], with Hi(0) = 1 and 
H 2 (0) = 0. Set 




(122) 



Xe(t,X,v) 
Xe(t,X,Z) 



^[e-H^H^E-^x - vt\] 



(123) 
(124) 




for (t,x) G M 2 . Then, for e > small enough, 




(125) 
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We apply Green's formula to the r.h.s. of (11251) . and let e — > and Z grow 
to 3V Using approximate limit property at (9, z) and (o) of Lemma 1X31 one 
obtains 

limsu P £- 1 (-/0[(M + £) x (0,1)] < V Tr[p(e,z-),p(e,z+)] 

Letting next (<5, a) — > (0, 0), we have 

A(0) < dd(9) + a(p + ,v(p-,p + )) + 7i(p~,p + ) (126) 

+ ]T tt[p(M-),p(M+)] 
zey e \{v} 

By (o)-(ii) of Lemma ELS we have A(0) < 0, and A(0) = implies ffTOTj) . 
ffT02D- flT03|) . and y[p(0, .)]\M = 0- Now, instead of (jI2BI) . we write 



-^-[(0,0 + e)x(O,l)] < - y Xe(t-9,x-y)dp(t,x) (127) 

- + x ((0,l)\(j/-We,y + We))] 

It follows that 

A(0) < <9d(0) + a(pV(p-,p + )) + 7r(p-,p+) 

- lim sup e~ V" [(0, + e) x ((0, l)\{y - We, y + We))) 

Thus A(6) =0 implies the following, which will be used in the last part of 
the proof: 

\mxe- l p-[(e,e + e) x ((0, l)\(y - We,y + We))] = (128) 
Third case, y = G y[p{9, .)}, Suppose first {9, 0) ^ J. By VMO property, 

/d+e t-We 
/ - Pol dxdu — (129) 

Jo 

We apply (1591 with ip = ips, £o = f'(po) > 0. Proceeding as in the first case 
one shows that A(9) is bounded above by (1931) . 
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Suppose next that (0,0) G J. Then (see Remark 13.1 j) (n t = 0,n x = 1) 
is the normal to J at (9, 0), and 

/9+e re 
/ |p(u,ar) - p(#, + ) I dxdw = (130) 
Jo 

If p(9,0 + ) < p*, we proceed as above, with £ = f'(p(0,0 + )) > 0, to show 
that A(6) is bounded above by (|93"j) . Finally, assume p(6>, + ) > p*. We then 
apply dHDD with y. = 0, thus d°(.) = and 9tZ(.) is given by ((HID. We bound 
the measure term as in (I125p . Since 9 G T, we arrive at 

A(0) < + ]T 7r(p(0,z-),p(M+)) (131) 

By (o), (ii) and (iv) of Lemma 13.3) we have A (9) < 0, and A(9) = implies 
(PSD, P5D and y[p(#, .)] n (0, 1) = 0. 

Fourth case, y = 1 G 3^[p(#, •)]• This case is similar to the third one, 
here we have to use (1891) with ■?/> = V'-S) or ASP)) with y. = 1. We arrive at 
ffTBTD . and A(0) = implies ([102]), ((Kg]) and y[p(fl, •)] n (°> x ) = 0- 

We now prove that, if ([§§!) is an equality, p~ vanishes on Q (which is open 
because of Lemma I'd. 51) . It follows from Lemma I'd. 51 that y. and y defined by 

y e := sup y[p(9,.)], y. := inf y\p(9, .)] (132) 



are respectively u.s.c. and l.s.c. By Proposition 13.1) they are continuous at 
every 9 G T, where we have y~g = y^. Applying ( )99|) on (u, t) with s < u < t 
shows that inequality (1991) is reversed on (s, w) and thus an equality, hence an 
equality any subinterval of (s,t). Let y, a Lipschitz-continuous [0, 1]- valued 
path on [a, b] C (s,t) such that y a G 3^[p(a, •)] and y& G 3^[p(&, •)]• Then, by 



(T((a,b) x (0,1)) < C(b-a) (133) 

where C is a constant depending only on / and h. Given rj > 0, we define a 
nonnegative measure M v < /i" C A on (0, +oo) by 

M V (A) = p- ({ (t, x) G A x (0, 1) : t G A, a; £ - r?, ft + 77] } 

for every Borel set A C (0, +00). Let 9 eT. Since and ?/ are continuous 
at 0, (PBBD implies 

lime- 1 M II ((0,0 + e)) = O (134) 
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Together with fl 133j) this implies M v = 0. Letting rj — > we obtain M = 0, 
i.e. pT vanishes on Q! := {(0,y) G (0, +oo) x (0, 1) : y £ [y g ,ye]}- Since 
fl\fl' is contained in U x (0, 1), where U is the negligible set of t G (— oo, 0) 
such that y < y t , using f)133p . we obtain /x~(f2) =0. □ 



Proof of proposition \3.2i Assume fllOOp is an equality. Then applying 
on (s, 0) shows that S[p(s, .)] > J(_oo,s) [p(-> ■)] which must be an equality in 
view of (DUD. Thus 

%(t2,.)]-%(*i>-)] = %,fe)b(- ) -)] (135) 

for every h < t 2 < 0. Let U := {t G (-oo,0) : .)] fl (0, 1) ^ 0} and 

*_ := inf U. We set t_ = if W = 0. By (II03J) - (jI05|) . the following occurs for 
a.e. t < t_: y[p(t, .)] = {1} if /(» < /(p r ), y[p(t, .)] = {0} if f( Pl ) > /(p r ), 
y[p(t, .)] C {0, 1} if f(pi) = f(p r )- By Lemma [3751 in the first (resp. second, 
resp. third) case, for all t < i_, 1 (resp. 0, resp. at least one element of 

{0,1}) lies in ;%(*,■)]■ 

First case. f(pi) ^ f(pr)- F° r notational simplicity, we treat the case 
fi.Pi) < f(Pi), which implies p x < p*. Then v := v(tp(pi),(p(p r )) < (cf. 
(|56|) ). Translation to the case f(pi) > f(p r ) is straightforward and left to 
the reader. By (jI35]) and ffTU^ - ffTptj) . for a.e. t < i_, we have p(t, l~) = p h 
p(t,0+) G £~(p ; ) and y[p(t,.)} = {!}. We apply §U§ on (-oo,t_] with 
= 1. By (ii), (iii) of Lemma [3731 we have d{t) = 0. Thus equality in f[5^1) 
implies p~[(— oo, £_) x (0, l)\ri)]=0, where Ti = (— oo, 0) x {1}. Hence p(., .) 
is an entropy solution to (T5]) on (— i_,+oo) x (0,1) and satisfies (1631 (here 
with jjoo = and cr = +). 

Now choose an arbitrary s G U and y G 7y[p(s, .)] fl (0, 1) and consider 
the upper forward solution y, of (11081) issued from (s,y). Let s' := inf{t G 
(s, 0) : ?/i G {0, 1}}, where the infimum is taken to be if the set is empty. 
We first argue that 

y t G y[p(t, .)] for every t G (s, s') (136) 

Indeed, assume there exists si G (s, s') such that y Sl G' y[p(si, .)). Let 
s := sup{£ G [s, si) : y t G [V[p(£, ■)]} ■ Th en s G [s,S i) , y So G y[p(s , •)], and 
for t G (s , si), (t, y t ) lies in the set fTTOTj) . By ffTUoD . ffl09l and Lemma [3761 
p(£, = p(t, yf) for a.e. t G (s , Si). We now apply flSUj) to p(., .) and on 
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[so, Si]. (i),(ii) °f Lemma [3~31 imply d(t) < 0, hence 

%0i,-)] - S\p(s ,.)] < S[p(si,.),y si \- S[p(s Q ,.),y S0 ] 

which contradicts (11351) . (1 10 II) and f[5"6l) imply y t = v(<p(p r ),(p(pi)) =: v a.e. 
on (s, s'). Hence y t = y s + v(t — s) for t G (s, s'). Since s was arbitrary, 
and (by Proposition 13. ip y[p(9, .)] is reduced to a single point for a.e. 9, the 
above implies i_ > — oo and existence of t + G (t-,0] such that y[p(t, .)] D 
{1 + «(*-*_)} for t G [*_,*+], and l+v(t + -t_) = if t+ < 0. By (jIDBjl and 
(HH3D, /i~|>(.,.)] vanishes on {(t,x) G (*_,*+) x (0, 1) : x ^ 1 + f (t - *_)}. 
Thus p(.,.) is an entropy solution to (j2|) on G (—£+,—£_) x (0,1) : 

x ^ «(t + £_)}. ffTOT]) and ([102|) ffT03|) of Proposition O yield ([ST]) and ® 
for p(.,.) on (-*+,-*_). 

We will now prove t + = 0, which will establish the proposition with y = vt_. 
The corresponding 9 y in Theorem 12.31 is — 1_. Assume £+ < 0, then by what 
precedes we have G y[p{t + , .)}. Define ti := inf{t G (t+,0) : y[p(t, .)] n 
(0, 1) 7^ 0}. Set t\ = if the preceding set is empty. By (I106p . p~[p(-, •)] van- 
ishes on D\ := (t + ,t 1 _) x (0, 1). Proceeding as above, if t l _ < 0, one can show 
existence of t\ G (*i,0] such that y[p(t, .)] D {l+v{t-ti)} for t G (*!,*+), 
and p~[p(.,.)] vanishes on _D 2 : = {(^j^) : x 7^ 1 + f(t — t 1 :)}. 

Let -D = £>i U D 2 . Hence, p(., .) is an entropy solution to fl2]) on the open 
set D := D\ U -D2, where .D, := $(-D«) is the space-time reversal (cf. (12 Op) 
of A- Besides, by f lT02|) - fri05|) . p satisfies (JSU on (-t%,-iL) and ([63]) on 
(— tl_, —t + ), with = and a = +. For e > 0, set 



where po < P*, an d define S" 6 as S in (|4T|) (148]) . but with y G [e, 1] and 
integral over [e, 1]. It is straightforward to translate ( 180]) to S E , a dynamical 
functional J e with left boundary e, and a [e, l]-valued path. Using (186 p . one 
can see that lim^o J £ [p £ (., .)] = J[p(.,.)]. Let (yf : t G [t+,0]) denote the 
upper forward solution of 



on [t+,0] issued from (t+,e). Since p(., .) is en entropy solution on Z), p(., .) 
has locally bounded space variation on D. In particular, if t G (t+,ii.) and 



p e (t,x) := p l(_oo, £ )(x) +p(t,x)l( £i+00 )(x) 



(137) 



y t = /'(p £ (*,m)) 



(138) 
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e < 1 + min(0,f(t — t^)), the limit p(t,e + ) exists in classical sense. Thus, 
(11371) and the Filippov sense of (11381) imply 

y t >e, p £ {t,e + ) > p* a.e. on {t : y t = e} (139) 

Note that, with the upper forward solution of y. of (11081) issued from (t + , 0), 
we could not deduce the analogous property p(t, + ) > p* on {y t = 0}, 
because the latter limit is only known in the sense (jSBJ. This is why we 
introduced th e-boundary. Set T £ = mi{t > t + : yf G B}, where B := 
[t+,tL) x {1} U {(t, l + v(t- £)) : t G [*!,*+]}. Since y £ has a Lipschitz 
constant ||/'||oo independent of e, T := liminf{T £ : e > 0} > 0. Let 
8 G (t+, T), T £ := {t G [t + , s] : 1 > y\ > e}, and T Q e := {t G [t+, s] : ^ £ = e}. 
By Lemma ESI P e {t,yf-) = P £ (t,y^ + ) a.e. on T £ . Applying (I5U|) to y E , using 
(I139p . (I103P and (i), (ii), (iv) of Lemma |3~3| we obtain 

S £ [p £ (s,.)]-S £ [p £ (t + ,.),s]< Jf t+ , r ,][p £ (.,.)] 

Since (y £ :,e > 0) has uniform Lipschitz constant it is relatively compact in 
C°([t + ,T}). Let y_ be a limit point, hence y t+ = 0. Since G y[p(t + ,0)], 
letting e — > yields 

%(«,•)> yj - S[p{t + ,.)} < J [t+)S] [p(.,.)] 

for all s G [£+,T]. This and (1 135|) imply y s G y\p(s, .)]. Since !V[p(s, .)] is a 
singleton for a.e. s and S is contained in the set (I98p . we obtain a contra- 
diction with the continuity of y_. 

Second case. f(pi) = f(p r )- By (I106p . p~[p(-, .)] vanishes on (— oo. t_) x (0. 1). 
For y G {0,1}, Let t y := inf{t G (-oo,t_) : y G y[p{t, .)]}. Set ^ = t_ if the 
corresponding set is empty. For at least one y G {0, 1} we have t y = — oo. 
We assume e.g. y = and leave the similar case y = 1 to the reader. Let 
r < £_ such that G !V[p(t, .)]. We proceed as above, now defining p £ , S e 
and J e with boundaries both at e and 1 — e. (I137P becomes 

p e (t, x) := p l(_ OOi6 )(ar) + p(t, x)l (£; i_ £) (x) + pil ( i i+oo) (x) (140) 

with p < p* < p\. Let y e be the upper forward solution of (I138P on [r, i_) 
issued from (r, e). Now ( 11 39[) becomes 

y| G [e, l-e], p £ (t,e + ) > p* on {yt = e}, p e (t, (l-e)~) < p* on {yf = 1-e} 

(141) 
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Let 8 G [r, T £ := {t G [r, s] : G (e, 1 - e)}, % £ := {t G [r, s] : y* = e}, 
7^ £ := {t G [r, s] : y\ = 1 — e}. Applying Lemma I3"l)l on T e and (150]) on [r, s], 
using (I14ip and (i)-(iv) of Lemma [3 .3[ we obtain 

S s [p%t-.),yL]-S £ [p £ (T,.),e}<J*[p*(.,.)} 

Let y, be a limit point of y e in C°([r, £_]) as e — > 0, then 

%0, 0>yJ -SKO] < j m[p(-,-)] 

In view of f)135p . this must be an equality with y s G y[p(s, .)]. Since y[p(s, .)] 
is a singleton for a.e. s and y, is continuous, we must have y s = G 3^[p(s, .)] 
for all s G [r, £_]. Since r is arbitrarily small, (I103p -( TT05|) imply that on 
(— 1_, +oo), p(., .) is an entropy solution that satisfies ( )63|) . here with ?/oo = 1 
and cr = — . 

If U = 0, i.e. i_ = 0, the above agument establishes the proposition with 
y — 1 — y. Assume W ^ I, L < 0. Since now we have v — 0, the same 
arguments as in the first case show that in fact U = (— oo, 0), £_ = — oo, and 
there exists y G (0, 1) such that p~[p(., .)] vanishes on (— oo, 0) x [(0, l)\{y}], 
p(t,0+) G £~(p*) and p(t, 1~) G £ + (p r ) for a.e. t G (-oo,0). Thus p(.,.) is 
an entropy solution on (0, +oo) x [(0, l)\{y '■— 1 — ?/}] that satisfies boundary 
conditions (ISTj) . This establishes the proposition with y = 1 — y. In both 
cases t_ = and i_ = — oo, the corresponding 9^ in Theorem 12.31 is — 1_. □ 

Proof of proposition \3.3{ Define p + (.,.), resp. p~(., .) on (0,0^) x (— oo, 1), 
resp. (0, 9 y ) x (0, +oo) by 

p + (t, .) := <p(Pr)l(~oo,y+^) +p(t,.)l (§+^,1) (142) 

p~{t,.) := ^(pj)l ( j( +trt)+0O ) +p(t, .)l( ,y+O (143) 

Since p ± has no jump across the line x = y + vt, is constant on one side, and 
p is an entropy solution to (J2J) outside this line, p ± is an entropy solution on 
its domain. Besides we have 

<P(Pr)l(-oo,y) + P(0, 1 - .)%!) (144) 

+ p(0,.l-)l (0>fi) (145) 
£~(p*) (146) 
£+(p r ) (147) 



P + (O,0 = 

p-(o,.) = 

p + (t,l~) G 

p-(t,0+) G 
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where (I146j) - (jl47p follow from (I6T]) . Thus p^ is uniquely determined as the 
entropy solution of an initial-boundary problem for ([2]). This uniquely de- 
termines p(., .) in (0, 6^) x (0, 1). Next, since pb G £ a (pb) for any o G {+, — }, 
f )63|) implies that on (6^, +oo) x (0, 1), p(., .) is the unique entropy solution to 
02]) with Cauchy datum p(6^, .) (previously determined) and boundary con- 
ditions p{t, 0+) e £ + { Pb ), p{t, 1-) e £~(p 6 ). 

Statement (0) of Theorem [273] follows because, assuming existence of /?(., .) : = 
M^[p(.)] for y 1 7^ y 2 , : if both y l G (0, 1), both p 1 have a single antishock 
along two two different curves. If exactly one y % G (0, 1), or if they both lie 
in {0, 1} and f(pi) ^ f(p r ), one p 1 has an antishock and the other does not. 
If both y % G {0, 1} and f(pi) = f(p r ), asumme e.g. y 1 = 0, y 2 = 1. Then 
p l (t,0 + ) = pi a.e. while p 2 (t,0 + ) G £ + {p r ) a-e. Since p\ < p* < p r , we have 
(see d38D (ESD) Pi £ + (p r ), which implies pV p 2 . 

Let us now prove that p(., .) given by f )57]) lies in 1Z T [p(.)} for some T > 0. 
If 9^ < +oo, this follows from Theorem 12.11 and pb ^ p* . Assume now 
f(Pl) = f{Pr), 9 y = +oo, y G (0,1). Since tp(pi) = p r G £~{p r ) and 
<p{Pr) — Pi G £ + (pi), by (I6T]) and (162]) . the restriction of p to (0, +oo) x (0, y) is 
an entropy solution to ([2]) with BLN boundary conditions p(t, + ) G £ + (p r ), 
p(t,y~) G £~(p r ). Similarly, the restriction of p to (0, +oo) x (y, 1) is an en- 
tropy solution to §2§ with boundary conditions p(t,y + ) G £ + (pi), p(t, l - ) G 
£~{pi). The result follows from Theorem 12.11 Pi P* > Pr P* ■ m all cases, 
we have lim^-oo .) = p 1_ * 00 (.), with p y (.) given by f l4~T]) . 

We next prove that p(.,.) given by fl57]) achieves equality in f l55]) . Set 
Zt = 1 — y-f Let — oo < T < 0. By (|80p applied to z, Lemma 13.31 and 
conditions (2), (a)-(c) of Theorem 12.31 

S\p(0,.),zo] -S[p(T,.),z t ] = J (i , )[p(.,.)] (148) 

Since p(0, .) = p(.) and 2 = y G ^[p(.)], f TT48]) and ([99]) on (T, 0) yield 
S[p(T, .), z T ] = S[p(T, .)], i.e. z t G y[p(T, .)], or equivalently y_ T G ^[p(T, .)]. 
For T small enough such that p(., .) G 7?/r[p(.)], (I148j) yields equality in ( 151?]) 
and fl55]) . 

We eventually prove that D^[p(s, .)] = {y s } for all s > 0. Assume the contrary 
and let i/ 1 7^ y 2 for some s > 0. Up to this point we have already established 
uniqueness of M^[p(.)] and statement (0) of Theorem 12.31 By anticipation 
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we use the existence of [/?(.)] established independently in the next sub- 
section. We set j?(t, .) = M^[p(s, .)](t - s, .) for t > s and p'(t, .) = p(t, .) 
for < t < s. Statement (0) implies p 1 ^ p 2 . But the uniqueness statement 
of Theorem 12.31 yields /?(., .) = p(., .) for each i G {1,2}. □ 

3.3 S > V and existence of a minimizer 

We construct, for each y G 3^[p(.)], a p(., .) satisfying conditions (2), (a)-(c) 
of Theorem 12.31 

Construction on the time interval (0, 9 y ). We proceed inductively by 
constructing p(., .) on intervals [t n ,t n+ i\. Let t = 0, z n = y tn = y + vt n , 



t n+ i -t n = min — — -± (149) 
\V — v V + v J 

where V = sup \ f'\ > \v\ by strict concavity of /. Then lim^oo^ = 6^. Set 
p(0, .) = p(.) and define p(., .) on (t n , t n+1 ) x (0, 1) by 

where p+ and p~ are the entropy solutions to (j2J), respectively on (— oo, 1) 
and (0, +00), with initial/boundary conditions 

Pn(M = <p(p r )l(-oo,^) +p(t n ,-) 1 ( 2 „,l) (151) 
Pn(O>0 = <P(P/)l(^,+oo) +p(t„,.)l( ,^) (152) 

p+(t,l-) G S-( Pl ) (153) 
p-(t,0 + ) G £+(p r ) (154) 

( ISTj) follows from fl 1 5 3 [) — fl 1 54 j) . We prove by induction on n that fl62|) holds 
a.e. on (0, t n ) and that 

z n e y[p{t n , .)] (155) 

For n = 0, (Il55p follows from y G 3^[p(.)]. Assume the statements hold for 
n. We first prove that ( )62|) holds up to t n+ x i-e. (in view of ( I150p ) 

p+(t - t„, (z„ + u(t - t n )) + ) = <P(Pr) (156) 
p~(t-t n , (z n + v(t-t n ))~) = ip(pi) (157) 
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for a.e. t G (t n ,t n+1 ). Let p+ + (.,.) and p n (., .) be the entropy solutions to 
@ on R with initial data 

p+ + (0,.) := p+(0,.)l(-oo,i)+Pcl(i,+oo) (158) 

Pn + (0,-) := /5^(.)l(0,+oo)+Pcl(-oo,0) (159) 

with p c given by flU}. By (fT55]) -(fT5D ]) . fTTBT]) - (jI52]) and (|I551» . 

p+ + (0,x)cte < p c (z-z n ), \/z>z n (160) 

p~~(0,:r)Gfe > p c (z n -z), \/z < z n (161) 



Given (TT58D fTT59|) . fTT60|) ( TT6T]) and (SSJ, we apply Lemma O below to p++ 
with r = <p(p r ) < P = Pc, and p~~, with r = <p(pz) > p = p c . Note that strict 
convexity of / implies 

v((p(p r ),p c ) > v((p(p r ),(f(pi)) =v 

v (r, p) = v((p(pi),p c ) < v(<p(p r ), (p(pi)) = v 

Thus Lemma I3T71 yields 

p + n + (t,(z n + v(t-t n )) + ) = ip(p r ) (162) 
p~~{t, (z n + v(t-t n )~) = ip{pi) (163) 

By finite propagation property (see e.g. Proposition (2.3.6) of [S]) for (T5]), 
we have p+(t, .) = p++(t,.) on (-oo, 1 -V{t-t n )) and p"(t, .) =p-~(t,.) on 
(F(t-i n ), +oo). With ffT62D- (TT63|) and ffT49D . this implies (TT56D on (i ft , t n+ i). 
Now we prove (11551) for n + 1. For t G [— t n +i, —t n ], set y 4 = 1 — z n — v(—t — 
t n ) G (0, 1), so 1 - z n+1 = y_ tn+1 . Define p(., .) by §7]), so (JISSD-dlSZD and 
(rr53l)- (IT5i|l imply ffTDTjl -fTTUSl) for a.e. 6 G (-t„+i, -*„)• Applying (JHOD with 
this y., using (o)-(ii) of Lemma [3TB1 we obtain 

S\p(-t n ,.),y- tn ] - S\p(-t n+ i,.),y- tn+1 \ = J(-t n+u -t n )[p(-,-)} (164) 

Then (199]) and (I155P imply (I155P for n + 1. Besides, summing (I164p over n 
and letting n — > oo in (11551) . we have 

%(.)]- %(-0~V)] = J ( _^ 0] [p(.,.)] (165) 
yoo e (166) 
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Construction on the time interval [9 y ,+oo). If 6^ < +00, which occurs 
unless f(pi) = f(p r ) and y G (0, 1), we define p on (6^, +00) as the unique 
entropy solution to ([2]) on (0, 1) with initial datum p(6^, .) and BLN boundary 
conditions 

pit, 0+) G £ + ( Pb ), pit, 1") g £-( Pb ) (167) 

with pb defined in (c) of Theorem 12. 31 The second boundary condition in 
(|63|) is thus satisfied. Let us show that the first one holds, i.e. 

p(t,r) = p b (168) 

with a and pb given in (c). We prove by induction that, for n G N and 
t n : = & + n /V, 

yoo ey[p(-t n ,-)} (169) 

and that this implies (11681) a.e. on (t n ,t n+1 ). For n = 0, (I169p is simply 
(11661) . We consider the case y^ = 0, as fjoo = 1 is symmetric. Here p b = pi, 
thus we must show 

p(t,0 + )= Pl (170) 

for a.e. t G (t„,t n+ i). Let p+(.,.) and p+ + (.,.) be the entropy solutions to 
([2]), respectively on (—00, 1) and R, with initial/boundary conditions 

Pn(°, = fan, -)l(0,l) + Pjl(-oo,0) (171) 
p+ + (0, .) = p(t n , .) 1(0,1) +Pll(-oo,0) +Pcl(l,+oo) (172) 

G (173) 



By (M), 



Note that 



p+ + (0,z)dz < p c x 



Voo = =S> f( P i) < f(p r ) ^Pi<p* (174) 

which, by fT4"9l . implies pi < p c . Thus Lemma I3T7I yields p+ + (t,x) = pi for ev- 
ery x < tv(pi, p c ). By (1491) and strict concavity of /, v (pi, p c ) > v(pi, y?(p r )) > 
0, where the last inequality follows from f{<p{p r )) = f(Pr) > f{pi)- Thus 

p++(i,0 + ) = p, (175) 

for a.e. t > 0. By finite propagation property, p + (t, .) and p ++ (t, .) coincide 
on (—00, 1 — V(t — t n )). Thus, by definition of t n , we also have 

plit^ + ) = p l e£ + i Pl ) (176) 
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for a.e. t G (0,t n+1 - t n ). By (fWj) . (fTTTj) . (fI73]) and ([175]), p(t n + .,.) and 
p+ are entropy solutions of a common initial-boundary problem for (j2J) on 
(0, t n+ i — t n ) x (0, 1). Thus they coincide on this domain. This in particular 
establishes (I170j) on (t n , t n +i)- We now apply ([50]) on [—t n +i, —t n ) with y t = 1. 
Using (IToTj) . ffTTOj) and f TTTl) . with (i), (ii), (iii) of Lemma E31 we obtain 
(HMD- Then ffToT)]) and (JHH]) imply ffToT)]) for n + 1. 

Lemma 3.7 

Let < r < p < K , and u(., .) 6e the entropy solution to |D) on R, with initial 
datum u (.) G L°°' K (M.) such that u (x) = r for x < and J* u (z)dz < px 

for x > (resp. u (z)dz > —rx for x < andu (x) = p for x > 0). Then 
u(t, x) = r for every x < tv(r, p) (resp. u(t, x) = p for every x > tv(r, p) ). 



Proof of lemma \J77\ We prove the first half of the statement, the other one 



being similar. We have 

u(t,x) = -d x U(t,x) (177) 
where U (t, x) denotes the viscosity solution to the Hamilton- Jacobi equation 



X 



d t U = f(-d x U), U (x) = -I u (z)dz (178) 
given by Hopf 's formula, 



o 



U (t, x) = inf 



u (y) + tf* 



x-y 



t 

where the convex conjugate of the concave function / is given by 



;i79) 



/» : = sup [f( P ) - vp] = mr 1 ^)] - vur\v) (iso) 

pe[o,K] 

and satisfies 

cry = -c/T 1 (isi) 

Set 

'x - y 



H(y) = U (y)+tf 
For y < we have, by assumption on w (.) 



t 



H( !J ) = -r !1 , If [^j 1 
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hence, by (11811) 

H'(y) = -r + (fT 1 (^) (182) 

Assuming x < tv (r, p), we have x — tf'(r) < 0, since concavity of / and r < p 
imply v(r,p) < f'(r). Thus H'(y) < for y < x — tf'(r) and H'(y) > for 
x — tf'{r) < y < 0, so that 

inf H(y) = H(x - tf'(r)) = -rx + tf(r) (183) 

y£(-oo,0) 

For y > we have, by assumption on u (.), 



H(y)>G(y):=-py + tr( 



* ( x y 



t 

Two cases may arise: (a) If x — tf'(p) < 0, then we have 

inf H(y) > inf G(y) = G(0) = H(0) > inf H(y) (184) 

j/£(0,+oo) j/G(0,+oo) j/e(-oo,0) 

(b) if a; — tf'(p) > 0, then we have 

inf H(y)> inf G(y) = G(x-tf'(p)) = -px+tf(p) > -rx+tf(r) = inf H(y) 

ye(0,+oo) ye(0,+oo) 2yG(-oo,0) 

(185) 

where the last inequality follows from x < tv(r,p). From (I183p and (j!84p - 
(11851) we obtain that, for every x < tv(r, p), U (t, x) = —rx + tf(r), and thus 
p(t,x) = r by (TT77D. □ 



4 The case p\ > p r 

4.1 Analysis of entropy production 

Lemma 4.1 Let p(., .) G £((— oo, 0) x (0, 1)). Let G (.,.), defined on (0, +oc)x 
M, 6e t/ie function defined in Theorem \2.J\ For (t,x) G (— oo,0) x (0, 1), set 
Ft(x) = ip[G(—t, 1 — x)]. Then, for every — oo < a < b < 0, 

S[p(b, .), F b (.)} - S[p(a, .), F a (.)] - J {a , b ) \p(; •)] = f d(t)dt - p-[(a, b) x (0, 1)] 

J a 

(186) 
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where d{t) = dd{t) + d°(t), with 

d°(t) := -d x F t (x)h"(F t (x)) (187) 

x {f(p(t, x))-fo <p(F t (x)) - f o <p[F t (x)](p(t, x) - <p[F t (x)])} 
dd(t) := g(p(t,0 + ),F t (0+))-g(p(t,r),F t (l-)) (188) 

Lemma 4.2 Assume p\ > F + > p* > F~ > p r . Then: 

1) For every p^ G [0, K], 

g{p + ,F + ) < i\p, Pl ) (189) 
-g{p-,F-) < f(p,p r ) (190) 

2) Equality occurs in M89\) iff.: either F + = p\ and p + G £~{p{); or, F + < p\ 
and p + G {F + ,(p(F + )} C £~(F + ). Equality occurs in U90\) iff.: either 
F~ = p r and p~ G S + (p r ); or, F~ > p r and p~ G {F~,(p(F~)} C £ + {F~) 



Proof of lemma 4-2 We will prove the statement for p + and F + , the other 
part being similar. The case F + = p\ follows from ii) of Lemma Assume 
now F + < pi. If p + > p* , then by Lemma [3.41 and (1351) . 

g(p + ,F+)<0 = i\p + , Pl ) 

with equality iff. p + = F + . We assume now p + < p*, and distinguish two 
cases, i) (p(pi) < p + < p*. Then by fl35l and Lemma [3.41 

g{p + 1 F + ) -i l (p + , Pl ) = g^(p + ),F + ) < 

with equality iff. y?(p + ) = F + . ii) p + < ip(pi), then by (|35|) . i\p + ,pi) = 
g(p + ,pi). To study g(p + ,F + ) — g(p + ,pi) we note that, for p* < r < tp(p + ), 

± g (p+,r) = h"(r)[f(r)-f(p+)]>0 

Thus, since p* < F + < pi < f(p + ), we have g(p + ,F + ) < g(p + ,pi) with 
equality iff. F + = pi, in which case p + G £~(pi). □ 

Proof of lemma \4-l\ We write 

0,^(0] = #(*)- j(*) + «(t) 
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where 

H(t) := [ h[p(t,x))dx, J(t) := f [p(t,x)h'(F t (x))]dx, «(t) := / k[F t (x)]dx 
Jo Jo Jo 

By the Gauss-Green formula for divergence-measure fields (|15j). 

H(b) - H{a) = f [g(p{t, 0+)) - g(p(t, 1"))] dt + p[(a, b) x (0, 1)] (191) 

J a 

Since (5(0,.) is nonincreasing, for t > 0, G is a solution of class C 1 con- 
structed by standard characteristics (see (11921) — (1193)) below). Thus by (T5]) 
and the generalized product differentiation rule for the divergence-measure 
field (p(t,x), f(p(t,x)) ([151 Theorem 3.1]), 

d t [p(t,x)h'(F t (x))}+d x [f(p(t,x))h'(F t (x))} =p(t,x)d t h'(F t (x))+f(p(t,x))d x h'(F t 

Thus h'(F t (x))(p(t, x), f(p(t, x)) is a divergence-measure field, and the Green's 
formula yields 

J(b)-J(a) = f [f(p(t,0 + )h>(F t (0+))-f( P (t,l-)h\F t (l~))}dt 



a 



b /•! 



+ / / [p(t,x)d t h'(F t (x)) + f(p(t,x))d x h'(F t (x))}dxdt 



a JO 

Since G(.,.) satisfies we have a t F t (x) = -f'{ip(F t {x))d x F t (x). By (g2 
d t k(F t ) = ip(Ft)h"(Ft)dtFt, thus the result follows by a simple computation, 
noting that (since / o <p = f) 

1 f(p(F t (x))ti'(F t (x))d x F t (x)dx = [h'(F t (l-))f(F t (l-))-g(F t (l-))} 

- [h'(F t (0 + ))f(F t (0 + )) — g(F t (0 + ))} 

□ 



4.2 S < V and uniqueness of a minimizer 

We need the following properties of (5(., .) defined in Theorem 12.41 Property 
1) below implies ( 1671) in Theorem 12.41 
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Lemma 4.3 

0) For all t > 0, G(t, .) is a nonincreasing function. 

1) For all t > and x G (0,1): a) (p(pi) < G(t,x) < <p(p r ); b) pi < p* 
implies G(t,l~) = <p(pz), pi > p* implies G(t,l~) < p* ; c) p r > p* implies 
G(t 1 + ) = if(pr), Pr < P* implies G(t,0 + ) > p* . 

2) For every x G [0, l],G(t, x) —> f(p s ) as t — > oo ; where p s is the uni- 
form density defined in ( f^Op . Outside the MC phase (i.e. the third case in 
fijty), there exists T > such that G(t, x) = p s for allt>T and x G [0, 1]. 

Proof of lemma \4^\ Since G(0, .) is nonincreasing and / concave, G(.,.) 
is shock-free and can be constructed by standard characteristics. For every 
(t, x) G (0, +oo) x M, there is a unique y = y(x, t) such that 

y + */'(G(0, y-)) <x<y + tf(G(0, y + )) (192) 

and we have 

G{t ) x) = U')- l {^) (193) 

All statements are simple consequences of fl 1 9 2h — (11 93 j) . □ 

Let p(.,.) G n[p{.)]. We apply f lT86D on (T,0). Then d°{t) < 0, because 
h is convex, F t is nonincreasing, and / is concave and dd(t) < by ([189])- 
(I190p . ii) of Lemma I3.3l and 1), b)-c) of Lemma 14. 31 (note that if, for instance, 
Pi < p* , (H89p does not follow from Lemma 14.21 but from ii) of Lemma |3.3[ 
because then we have F + = F(t, + ) = pi by 1), b) of Lemma |4T3"1 p r > p* is 
treated similarly). Since F is the maximizing F for p(.) in floT]) . we have 

%(.)] < S[p(T, .), F T (.)] + J (m [p(., .)] < S[p T (.)} + J (Ti0 )[p(., •)] (194) 

Sending T — > — oo, since S'[p s (.)] = 0, we obtain 5*[p(.)] < J(_ OOj0 )[p(., •)], 
hence 5*[p(.)] < V[p(.)]. Let us now assume that p(.,.) G lZ[p(.,.)} achieves 
equality in fl5"4"|) . Then in (I186p we must have d°(t) = dd(t) = a.e. on 
(—oo,0). Thus 

S[p(t,.),F t ] = J { ^ t) [p(.,.)}>S[p(t,.)} 
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Hence F t = F p{u)l that is ([70]). Define p(., .) by (j57j). d°(t) = 0d(t) = on 
(—00, 0) has the following implications: 



1) Inequalities H189|) — fllQOH must be equalities for a.e. t < 0, with p 1 * 1 = 
p(t, O 1 * 1 ). Then 2) of Lemma 14.21 ii) of Lemma 13.31 and 1), b)-c) of Lemma 
EH imply p(t,0+) G £ + [G(t,0 + )] and p(t, 1~) G £~[G(t, 1")] for every t > 
(here again, if p^ < p* or p, r > p*, Lemma I3T31 is used instead of Lemma I4T2]) . 

2) p[(-oo,0) x (0,1)] = p+[(-oo,0) x (0,1)], hence p = p h [p(., .)] > 
on (-00, 0) x (0, 1), and by (J2U]), p A [p(., .)] < on (0, +00) x (0, 1). Since © 
is invariant by space-time reversal, p is also a weak solution. Thus, by [18], 
p(., .) is an entropy solution to ([2]) on (0, +00) x (0, 1). 

Since p(., .) G C°((-oo, 0), L^O, 1)), 1) and 2) imply that p(., .) is the en- 
tropy solution to 02]) with BLN boundary conditions ( 169]) and initial datum 
f ]68]) . Thus, if p G 7£[p(.)] achieves equality in (154"]) . it is necessarily (and 
uniquely) determined by the construction of Theorem 12.41 an d satisfies (|70|) . 

4.3 5 > V and existence of a minimizer 

Let us first recall the construction of the maximizing function F(.) in (150]) . 
Let R(.) be a continuous function on [0, 1]. and —00 < a < (3 < +00. We 
define the (a, /3)-truncated convex hull C a ^R of R(.) as the upper envelope 
of all convex functions S(.) on [0, 1] such that S < R and a < S"(0 + ) < 
>S"(1~) — P- A straightforward extension of a computation of [23] shows that 
the supremum in ( ]5T]) is achieved uniquely for F = F p (.) given by 

v[Fp(.)(x)} = -^C v(pi)Mpr) R(x) (195) 

where R(x) := J* p(y)dy. Alternatively, we have 

C a;f3 R(x)= sup [6x-R*(6)] 

fle[a,/3]nK 

where 

:= sup [0a:-i2(z)] 
xe[o,i] 

is the Legendre-Fenchel transform of R(.). C^oo^ooR = R** is the usual 
(untruncated) convex hull. One-sided derivatives (C Q , i ^i?)'(x ± ) G R are such 
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that [{C a ^R)\x ), (C at pR)' '(x + )} is the set of maximizers of the concave func- 
tion 9 i-V 9x - iT(0) on [a, (3}. Note that (C a ^ J R) , (x ± ) = {CR)'{x ± f a , where 
:= max[a, min(x, The above definitions are easily adapted to the 
truncated concave hull C a ^ for a > ft, by replacing suprema with infima and 
maximizers with minimizers. 

We shall need a Hopf-Lax type formula derived in [2H] for the entropy solu- 
tion to $5]) on (0, +oo) x (0,1) with Cauchy datum p Q (.) G L°°((0, 1)) and 
time-dependent boundary conditions 

p(t,0 + ) g £ + [ P i{t)l P (t,l-) G £-[p r (t)\ (196) 

a.e. in the sense of [3], where pi{.) and p r {.) are boundary data in L°°' K (0, +oo). 

Theorem 4.1 (JMfl) For x, y G [0, 1] and t > 0, let C(y, x, t) denote the set 
of differentiable paths 7 : [0,t] — > [0, 1] such that 7(0) = y, 7(1) = x and 
f'{K) < % < /' (0) a.e. Let R(0,x) = /* p{0,y)dy and define 

R(t,x) := sup \R(0,y)- l {0 }(7 s )/[min(pi(s), p*)]ds 

ye[0,l],y.eC(y,x,t) I JO 
t 

l{i}(7s)/[max(p r (s),p*)]rfs 
+ £l ( o,i)(7 s )/*a)^| (197) 

where f*(9) := mf pe [ 0t K] [0p — f(p)] is the (concave) Fenchel-Legendre trans- 
form of f . Then d x R(t,x) exists a.e. and is the entropy solution to (TJ|) with 
initial data p(0, .) and boundary data U96]) . Besides, the supremum in \19T\) 
can be restricted to the subset ofC(y,x,t) consisting of piecewise linear paths 
with constant slope between the boundaries. 

We now show that p(.,.) G lZ[p{.)} constructed in Theorem 12.41 achieves 
equality in (152]) and f )54|) . and that outside the MC phase it lies in TZt[p(-, •)] 
(thus achieving equality in fl53|) and f )55|) ). We interpret the restriction of 
G(., .) to (0, +00) x (0, 1) as the entropy solution to §2§ on (0, 1) with carefully 
chosen BLN boundary data and use Theorem 14. 11 Set R(0, x) := p(0, y)dy. 
By ([HS]), we have that G(0,x) = £§(0,x), where 

S{0, .) = C^'^RiO, .) (198) 
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Following flMD-flSSD, it is trivial that G(t, 0+) G £ + [G(t,0+)\ and G(t, 1") G 
l - )]. On the other hand, using 1), b)-c) of Lemma FO| one can 
see that if p r < p*, then we also have G(t, + ) G £ + [<p(G(t, + ))]. Likewise 
if p; > p* we also have G(t, 1~) G £~[<p(G(t, 1~))]. Choosing the suitable 
boundary condition in each case, T heorem 14.11 yields G(t, .) = d x S(t, .), where 

S(t,x):= sup \S(0,y)- l {0tl} (j s )f[G(s,j s )]ds + l(o,i)(7«)/*(%) ds 

2/e[Q,l],7.eC(j/,a;,t) I Jo JO 

(199) 

On the other hand, by definition, G(t, x) is the entropy solution to (j^J) on 
(0, +oo) x R with initial datum fl66|) . It is thus also given by the standard 
Lax-Hopf formula, hence for x G [0, 1], 

S(t,x) := sup \S (y)+ l(o,i){ls)f*(js)d \ 

ye[o,i],f.eC(y,x,t) I Jo J 

= sup !s (y) + tr(^-)\ (200) 

where 5(0, x) = L G(0,y)dy, and C(y,x,t) is the set of different iable paths 
7. : [0,t] -> R such that ^o = V, It = x and /'(if) < 7 S < /'(0) a.e. on 
[0, £]. (I200p has a unique minimizer y — y(x,t) given by (I192p . and satisfying 
( 11931) . Applying again Theorem 14. II to p(., .) defined in Theorem 12.41 we have 
p(t,x) = d x R(t,x), with 

R(t,x) := sup \R{0,y)- f l {0} (7 s )/[min(^[G( S ,0 + )],p*)]rf S 

ye[Q,l]n.eC(y,x,t) I JO 

- / l{i}(7s)/[max(<p[G(s, l~)],p*)]cte 
./o 

+ ^ l(o,i)(7.)r(7-)d«} (201) 

We are now ready to establish the following properties, from which the result 
of this subsection will follow easily. 

Lemma 4.4 For every t > 0, fiTffl holds, or equivalently, 

S{t, .) = C^^Rit, .) (202) 
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Corollary 4.1 For a.e. t < 0, equality holds in /I189\) - flfflfy , with p A 
p{t, 0+), p~ = p(t, F+ = F t (0+), F~ = F t {\-). 



Proof of corollary \4-l\ We treat (11891) . the argument being similar for ( I190p . 



Assume first that p\ < p*, then 1) of Lemma [4.31 implies ip(G(—t, 1~)) = p\. 
The result then follows from pit, + ) = p(—t, 1~) G £~[<p(G(— t, 1~))] and ii) 
of Lemma 13 .31 Assume now pi > p*. If G(—t, 1~) = (f(pi), we argue as above. 
If G(-t,l~) > cp(pi), Lemma S31 implies p(t,0+) = p(-t, 1~) > G(-t,l~). 
On the other hand, by 1), b) of Lemma 14.31 (p[G{— t, 1~)] > p*. This 
and p(t,0+) = p(-t,l~) G £-[<p(G(-t,l-))] imply (see (J3DJ) p(t,0+) G 
{(^(G(-t,l-)),G(-t,l-)} = {F t (0+),<^(F t (0 + ))}. The result then follows 
from 2) of lemma 14.21 □ 



Proof of lemma \4-4\ We check properties (l)-(3) below, which imply the 
result. 

1) S(t, .) is a concave function, y?(p r ) > d x S(t, + ) > d x S(t,l~) > (f(pi), 
and 

R(t,x) <S(t,x) (203) 

The first two properties follow from G(t, .) = d x G(t, .) and 0) of Lemma [4.31 
f l203|) follows from (JTjgHMJ, ^(0, .) < 5(0, .) (which is implied by (ELSE]) ), 
and the inequalities (cf. ( 133]) ) 

/[min(p[G( S , 0+)], p*)] > 0+)), f[msx(cp[G(s, 1")], p*)] > 1")) 

(2) ^(t, .) = .) on (0, l)\/t, where J 4 denotes the union of all (relatively) 
open subintervals of (0, 1) on which S(t, .) is affine. To prove this we write 
I t = (0, 1) D It, where I t denotes the union of all open subintervals of K. on 
which Sit, .) is affine. By (HH2D (HHSD ) , 

I t = {x eR: yix,t) el } (204) 

Assume now that x G (0, We first claim that y = y(x,t) G [0, l]\io- 

Indeed, suppose e.g. that y < 0, then f )66|) and H192|) — fll93|) imply that 
£ = y + tf'((p(pi)), hence p/ > p* , and G(t, .) = <p(p«) i n a neighborhood of 
x, which would imply x G It- A similar argument holds if we suppose y > 1. 
From ?/ G [0, l]\io an d (|198|) . it follows that 5(0, y) = -R(0, jy). The maximiz- 
ing path in (I200p is 7 S = y + sf'(p), with p = G(t,x) given by (j!93j) . This 
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path lies in C(y,x,t), thus the same y and 7. also produce the maximum in 
(11991) . Since 7. does not see the boundaries (except possibly at time 0), it 
produces the same value in (11991) and (12011) . Thus, in view of (12031) . y and 
7. must achieve the supremum also in (120 ip . and R(t,x) = S(t,x). 

(3) (a) If d x S(t,0+) < <p(p r ),thenS(t,0+) = R(t, 0+), and (b) Ud x S(t,l~) > 
(f(p r ), then S(t, 1~) = R(t, To prove this we set 

y := sup{x > : (5(0, x~) = <p(p r )} 
Vl := infix < 1: G(0, x + ) = <p( Pl )} 

Thus G(0,y -) = V { P r) > G(0,y +) and G(0, Vl -) > G(0, yi +) = V { Pl ). 

Let t = yo/f(<P(Pr))~ and h = (1 - Vi)/ 'f(<f(Pi)) + , witn tne convention 
0/0 = +00. Hence to = +00 iff. p r > p*, and t\ = +00 iff. pi < p*. One can 
see from (TT921) ffT93|) that 

t := inf{t > : G(t,0 + ) < <p(pr)} 
tx := inf{t>0: G{t, 1~) > cp{ Pl )} 

and, if to (resp. ti) is finite 

G(t ,0 + ) = cp(p r )>G{t,x), Vx> (205) 
G(t ,0 + ) = (pfa) < G{t,x), Vx> (206) 

In particular, we have nothing to check for x = and t < to, resp. x = 1 
and t < t\. If to = ti = +00 we are done. Otherwise, we assume e.g. 
to < t\. The above arguments imply that condition 3) is checked (since void) 
for t G (0,t ). Thus for such t, (I202p holds. Since p(.,.) and G(., .) lie in 
C°((0, +00), L^O, 1)), (JH2D holds also at t = t . We claim that 

£(t ,0 + ) = 5(t ,0 + ) (207) 

Indeed, assume i?(t ,0 + ) < S(t ,0 + ). Then there exists x > such that 
R(t ,x) < S(t ,x) for x G (0,x ). This implies S(t , .) is affine, and thus 
G(t , •) constant, on (0,x ), in contradiction with (12051) . Let t > t . We 
consider (I20ip with time origin at t and initial condition R(t ,.). Since 
t < +00 we have (by 1), b) of Lemma ETBl G(t, + ) > p*. Then the special 
path 7 S = for s G [t , t] in (12011) yields 

R(t,0 + ) > i?(t ,0 + ) - / f(G(s,0 + ))ds 
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On the other hand, S is the viscosity solution of the Hamilton Jacobi equation 

d t S(t, x) + f(d x S(t, x)) = (208) 

Integrating ( I208p on the space interval (0, +oo) and then on the time interval 
\po,t], and using G = d x S, we obtain 

3{t, + ) - S(t , + ) = - f f(G(s, + ))ds 

and thus (by f|2D7|> ) S(t,0 + ) < R(t,0+). With (F2DSD, this implies S(t,0 + ) = 
R(t,0 + ) for all t > t . Thus condition 3) holds on [to? ^i)- It follows that 
( 12021) holds on the time interval [to,ti). If t\ = +oo we are done. Otherwise, 
by continuity as above, f !202j) holds at t\. Then we proceed as above to show 
that (I207p holds at time ti, as well as 

R(t 1 ,l-) = S{t 1 ,l-) (209) 

and obtain R(t, 1~) = S(t, 1~) for all t > t±, whence condition 3) and (I202|) 
on this interval. □ 

We now conclude the proof of (I5"4"j) - (I55I) and Theorem 12.41 as follows. 

1) p(.,.) defined by fl57|) achieves equality in fl52l) and fl54l) . and equality in 
f )53|) and fl55|) outside the MC phase. Indeed, we apply (11861) on (T, 0) with 
T < 0. By flU}, At~[p(t, .)] = 0. Lemma S3] implies d°(t) = and Corollary 
SH] implies = 0. Let F(., .) be defined by (!7Dp. By Lemma S3 •) 
achieves supremum in ( 151]) for p(t, .). Thus ( I186P yields 

S[p(.)]-S[p(T,.)] = J im [p(;-)} (210) 
for every T > 0, hence SfpQ] > J(-oo,o) [p(-, •)]• 

ii) Outside the MC phase, p(., .) G TZr[p(-)] for some T > 0. Indeed, by 

2) of Lemma I4.3[ there exists r > such that, for t > r, p(., .) is the en- 
tropy solution to 02]) with Cauchy datum p(r, .) and boundary conditions 
p(t,0 + ) G ^ + (p s ), p(t, 1~) G £~(p a ). By Theorem [ZTJ since p s 7^ p*, p(.,.) 
relaxes in finite time to the steady state p s . 
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iii) In the MC phase, we still have p(., .) G lZ[p{.)}. Indeed, let e < min(p*, A— 
p*). By 2) of Lemma|131 we can find r > such that p* -e < <p(G(t, + )) < 
p* + e and p* - e < <p(G(t, 1~)) < p* + e for t > r. For t > r, let p+(t, .) be 
the entropy solution to (T5]) with Cauchy datum p + (r, .) = p(r, .) at time r, 
and boundary conditions p + (t, + ) G £ + (p* +e), p + (t, 1") G£ - (p* + £:). Sim- 
ilarly, define p~(t, .) with boundary data p* — e. By monotonicity of (J2J) with 
respect to boundary conditions, p~(t, .) < p(t, .) < .). By Theorem 12.11 
since p* ±e ^ p*, there exists r' > r such that /^(t, .) = p* ± e for all t > r' . 
Thus, for t > t', p* — e < p(t, .) < p* + e. This establishes pointwise (and 
L 1 ) convergence of p(t, .) to p s = p* as t — > +cxd, hence p(., .) G 72.[p(.)]. 
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A Variational expression for the dynamic ac- 
tion functional 

Recall from [4] that (F,G) is called an entropy- flux sampler for (J2]), where 
F(t, x, p) and G(t, x, p) are functions of class C 2 with compact support in 
(a, b) x [0,1] x [0, A], iff., for each (t,x), (F(t, x, .), G(t, x, .)) is an entropy- 
flux pair. We add boundary conditions ([3]) to this definition as follows. We 
say (F, G) is boundary compatible iff., for x G {0, 1}, (F(t, x, .), G(t, x, .)) is 
a boundary entropy- flux pair in the sense of [38J, i.e. 

F(t, x, p b (x)) = d p F{t, x, p b {x)) = (211) 

G(t, x, Pb (x)) = d p G(t, x, p b (x)) = (212) 

where the boundary data p b (.) is given here by p b (.) = pjl{o} + p r l{i}- The 
F-sampled entropy production of p(., .) G L°° ,A: ((a, b) x (0, 1)) is defined ([I]) 

as 

Pfa-,)-=- [ [ [{d t F)(t,x,p{t,x)) + {d x G){t,x,p(t,x)))dxdt (213) 

J a JO 
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Let /ibea given uniformly convex entropy. We denote by 7-\(0, 1) the set of 
boundary compatible entropy samplers F such that d p F(t,x, p) < h"(p). 

Theorem A.l For every p(., .) £ L°°' K ((a,b) x (0, 1)), 

7 (a , 6) [p(, •)] = sup{P F , p( . )0 : F £ P,(0, 1)} (214) 

Corollary A.l /or ever?/ c £ [0, +oo) 7 77 oo, c]) zs compact with respect 
to the local L 1 topology on L°°' K ((a,b) x (0, 1)). 



Proof of corollary \A.1[ For each F £ 7\(0, 1), p(., .) >-)■ Pp iP is continu- 
ous, thus 7( 0) 6) is lower-semicontinuous. We show that (7? b ))~ 1 ((— oo, c]) D 
(-^(a,6))~ 1 (( — oo, c]) is relatively compact. Consider functions 77 £ C 2 ([0,TT]) 
not necessarily assumed convex, that we call entropy by extension. The en- 
tropy flux (by extension) is still defined by (JTTj) . We may extend ( !T9|) and ( )2T|) 
to such functions (write f] = f]i — t]2, where rji are convex functions of class 
C 2 ). Let VL be a bounded open subset of (a, b) x (0, 1). Since h" is uniformly 
convex, by (I72"|) . .)](0) is bounded on (7° a ^) _1 ((— 00, c]). Then the 

compensated compactness argument following [44j Proposition 9.2.2] yields 
the desired result (notice that Tartar's equation (HI (9.15)] is stable by uni- 
form convergence of entropies, and thus holds for C 1 extended entropies 
including /). □ 

For the proof of Theorem IA.1| we need the 

Lemma A.l Let $(po) denote the set of entropy-flux pairs (i],q) such that 
v(Po) = v'(Po) = 0, q(p ) = q'(po) = (215) 

Then 

i\p, Pl ) = sup{g(p): (77,9) V"<h"} (216) 

i\p )Pr ) = sup{-g(p): (r/,g)£$(p r ), r)" < h"} (217) 

The suprema in \216\) - {211ty are respectively achieved by the following en- 
tropies 7] l pi p and rf pr p . 



rf pipi {u) := / (u-v)h"(v)l {0t+oo) [q v (p, Pl )]dv (218) 

J Pi 

Vp, Pr (u) ■= / (u-v)h"(v)l (0i+oo) [-q v (p,p r )}dv (219) 
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Proof of lemma DOl The result follows from f )271) applied to (i],q) G $(po)- 
To maximize q(p) and achieve i l (p,pi), one has to choose the entropy rj G 
$(p z ) so that 77" (v) = /i"(w)l (0 , +oo )[g t ,(p,Pi)] ) and similarly for i r (p,p r ). □ 

Proof of theorem \A.1[ 

Step one. We prove that the r.h.s. of (I214p if finite iff. p(.,.) G £((a,b) x 
(0,1)), and is then dominated by the l.h.s. Assume first that p(.,.) G 
£((a,b) x (0,1)). Then I( a> fy[p(., .)] is finite because t> m(v;dt,dx) is 
bounded. Using the generalized Green's formula ([IS]), one obtains the fol- 
lowing boundary extension of the formula established in [U Proposition 2.3]: 



Inequality > in (12141) then follows from F"(t, x, v) < h"{v), G(t, 0, .) G ® h {pi), 
G(t, 1, .) G &h{Pr) and Lemma [A. 11 Assume now that the r.h.s. of (12141) is 
finite. We first show that p(., .) must be a weak solution of ([2]). Assuming 
the contrary, there exists (p G C^((a,b) x (0,1)) such that Pp >p 7^ 0, with 
F(t,x,p) = Cip(t,x)p. Then the r.h.s. of (12141) dominates the supremum 
over entropy samplers CF(t,x, p), where C G M, which is +00. Now, fix a 
C 2 entropy r] such that rj" < h". Let F(t,x,p) = ip(t,x)r](p), where ip varies 
in C|?((a, b) x (0, 1)), < ip < 1. Then P F<p = p v [p(., -)](<p) remains bounded 
over ip, i.e. p^[p(., .)] G M((a, b) x (0, 1)). Since h" > c > 0, this is true for 
any C 2 entropy 77. Thus p(., .) G £((a, b) x (0, 1)). 

S'tep two. Assuming p(., .) G £((a,b) x (0,1)), we show that the r.h.s. of 
(12141) dominates the l.h.s. Let e > and a e : [0, 1] — > [0, 1] be a smooth cut- 
off function such that a e (x) = 1 for x G [0, e/2] and a e (x) = for x G [e, 1]. 
Set (3 e (x) = a £ {l — x) and 7 e = 1 — a £ — /3 e . We define a regularized (w.r.t. 
p) version of ff Pr P in (12 18[) by convolution 




(220) 



(221) 




(222) 
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where x is a standard convolution kernel. To make the above meaningful we 
naturally extend rf for p £ [0, K] by setting rf = vLk f° r P > rf = 
rf pifi for p < 0. Following f l2T8|) - fl2T9|) and (j2g])-(p|), and 7?^ p (V)'are 

ii-Lipschitz functions of u with ii = h'(K) — h'(0), and continuous (hence 
uniformly continuous) functions of {p,u). Thus convergences 

7£» »&», Cp(«) ( 223 ) 

are uniform with respect to (p, u) G [0,if] 2 . We also regularize boundary 
traces p(t, + ) and p(t, 1 _ ) by convolution, setting 

p E (t,0 + ) := J p(t-s,0 + ) X (s/e)ds (224) 

/f(M") := f p(t-8,l-) X (8/e)ds (225) 



We may view the mapping v h-> m(t>; c?t, dx) as a measure M(dt, dx, dv) = 
m(t>; (it, dx)dv on [0, if] x (a, 6) x (0, 1). Then M + (dt, dx, dv) = m + (t>; dt, dx)dv. 
Let ip e G Cf((a, b) x (e, 1 - e) x [0, if]), < ^ £ < 1, such that 

"K r f K f 

il) £ {t,x,v)m(v]dt,dx)dv > / / m + (v;dt,dx)dv — e 

'0 J(a,b)x(Q,l) JO J (a,b)x(e,l-e) 

(226) 

We define F £ (t,x, .) as the unique entropy (given by ( 12T|) ) r^(.) G $(po) such 
that 7/'(t>) = h"(v)i/) e (t,x,v), where po G [0,if] is arbitrarily chosen (but 
independent of (t, x)). Let 5 > 0. We apply (I220p to the entropy sampler 

■Su \ ._ / \ i t — a\ _ IS 



F*(t,x,p) := « £ W7 E ^JC j(i ,o + )W (227) 

+ 7e (a;) J P e (t,a;, i o) (229) 
We use (I226p . the triangle inequality 

\Vptp(p + ) ~ ^(P + , Pi) I < \v P ip(P + ) - Vp'tpiP) I 

+ \Vpip{p)-i l {p,Pl)\ 
+ \i l (p,Pl)-i l (p + ,Pi)\ 
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(note that r(p, p ; ) = T] m p(p)) with p = p 5 (t,0 + ), p + = p(t,0 + ), and the 



similar decomposition with i r (p ,p r ), p = p(t, 1 ) and p = p s (t, 1 ). From 
these we obtain that, in the limit lim^o lim^o, Pf s > s ,p{.,.) converges to the 
l.h.s. of ([211]). □ 

B Explicit minimizers for uniform p(.) when 

Pi < Pr 



In this and the next section, explicit minimizers are computed by interaction 
of Riemann waves, see e.g. [36] for details on such computations. Let p c be 
given by f l4"9]) . We consider p(.) = p, for some constant p 6 [0, K]. 

(1) Assume pi < p r and f(pi) < f(p r ) (thus pi < p*). Then: 

(a) If p < p h then y = 1, = , T v = l/f'( Pl ), p s (.) = p h 
Pi if x<tf'(p) 



p(t,x) = < 



(fTHx/t) if mm(tf( Pl ),l) < x <mm(tf'(p),l) (230) 



^ p if min(t/'(p)> 1) < x 

(b) If pi < p < p* and p < p c , then y = 1, ^ = 0, r» = l/v(p,pi), p s (.) = p h 

p{t, .) = pjl(o,Mp lW )) + * e [0, r y ] (231) 

with t>(p, pi) > 0. 

(c) If p* < p < p c , then y = 1, 0y = 0, r y will be defined below, p s (.) = p^, 

pi if < a; < tv(p,pi) 

p(t,x) = lp if tv(p, Pl ) <x< l + tf'(p) (232) 

k (f)- l ((x-l)/t) if l+f/'(p)<x<l 

fort<[t;(p,,p)-f(Pr 1 :=ii, 

p; if < x < x t 

U')-\(x-l)/t) if x t <x<l 
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p{t,x) 



(233) 



for ti < t < r y , where x t is defined for t > t x by 

x tl =t 1 v(p,p l ) = l + tj'ip), x t = v[p h (f')'\(x t -l)/t)] >v(p, Pl ) >0 

(234) 

and r y is the time at which x t = 1. 

(d) If p c < p < p r and f(pi) = f(p r ), then y = 0, 9 y = 0, r y = -l/v(p,p r ), 

Ps(-) = Pr, 

p(t, X) = pl(0,l+tv(p,p r )) + PA(l+tv(p,p r ),l) ( 2 35) 

with v(p, p r ) < 0. 

(e) If p r < p and f( Pl ) = f(p r ), then y = 0, 6 y = 0, r y = -l/f'(p r ), 

Ps(-) = Pr, 

p r if max(l + tf'(pr), 0) < x < 1 

p(t,a;) = i - l)/t) if max(l +tf'(p),0) < x < max(l + tf (p r ),0) 

p if x < max(l + tf'(p), 0) 

(236) 

(f) If p > p c and pi < p r < p*, then y = 0, 9 y = -l/v(ip(p r ),(p(pi)), r y is 
defined below, p s (.) = p h 

< {fY\{x-\)lt) if l + t/'(^))<^<l 

(f(p r ) if 1 + tv{ip{p r ),ip{pi)) < X < 1 +tf((p(p r )) 



p(t,x) = < 



¥>(Pi) 



if [1 + tv( P , ^{piW < X < [1 + tv((p( Pr ), ^{plW 

if < x < [l + tv(p,<p(pi))} + 



for t < d y , 



p(t,x) = < 



(237) 



Pi 
P 



if X < (t-0»)v(pi,<p(pr)) 

if t^(pi, (p(p r )) < X < 1 + tf'((f(p r )) 



( (f)-\(x-l)/t) if l + tf'(<p(p r ))<x<l 



(238) 
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for 9 y < t < h :-- 



v(f>l,<p(Pr)) ~ f'(<P(Pr)Y 

Pi if < x < x t 

(f)-\( Xt -l)/t) if x t <x<l 
for ti < t < r y , where x t is defined for t > t\ by 



p(t,x) = 



(239) 



X tl = l+ hf'faipr)), X t =V [ Ph (/') \{x t ~ l)/t)] > V{ PU (fi(p r )) > 

(240) 

and r y is the time at which x t — 1. 

(g) If p > p c and pi < p* < p r , then y = 0, = -l/v(<p(p r ), <p(p ; )), 

Ty = & + l/v(p l ,<p(p r )),p 8 (.) = Pl , 

(f(p r ) if 1 + tv(ip(p r ), (fi(pi)) < X < 1 

p(t,x) = ^ ¥>(p,) if [1 + to(p, <x< [l + tt;(v9(p r ),v?(Pi))] + 

^ p if 0<a;< [l + to(p )V 9(pO)] + 

for i < 0», 

P(t,x) = Pll(0MPlMPr))( X ) + <P{pr)Mtv{ Pl w{ Pr ),l){x) 

for §y <t< T y . 



(241) 
(242) 



(h) If p = p c < p*, and f(p t ) = f(p r ), then y G [0, 1] is arbitrary, 
r y is defined below, p s (.) = p|(.), 



= +oo, 



p(t,x) 



p c 


if 


min(l - y + tv(pi, p c ), 1) < x < 1 


Pi 


if 


1 - y < x < min(l - y + tu(pj, p c ), 1) 


Pr 


if 


1 - y + tv(p c , p r ) < x < 1 -y 


Pc 


if 


£f (Pc) < a; < 1 - y + tu(/0 c , Pr) 


(/rv/*) 


if 


< x < t/'(p c ) 



(243) 
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for < t < h : = (1 - y)(f'(p c ) ~ v(p c , p r )), 

Pc if min(l - y + tv(pi,p c ), 1) < x < 1 

pi if 1 ~y < x < min(l - y + tv(p h p c ),l) 

p r if x t < x < 1 — y 

{ {f')-\xlt) if 0<x<x t 



p(t,x) = < 



(244) 



for ti < t < t v , where r y = max(t 2 , y/v(pi, p c )), t 2 is the time at which 
x t = 0, and x t is defined for t > t\ by 



X h = tlf(p c ), X t =V [(f) \x t /t),Pr] < V(p c ,p r ) < 



(245) 



p(t,x) 



(i) If p = p c > p*, and f(pi) = f(p r ), then y G [0, 1] is arbitrary, 6 y = +00, 
r y is defined below, p s (.) = p y (.), 

{ {f')-\(x-l)/t) if l+tf'( Pc )<x<l 

Pc if 1 - V + tv(p h p c ) < x < 1 + tf'(p c ) 

pi if 1 ~y < x < 1 -y + tv(pi,p c ) 

p r if max(l — y + tv(p c , p r ), 0) < x < 1 — y 

p c if < x < max(l — y + tv(p r , p c ), 0) 

(246) 

for < t < h := y/(v(p h p c ) - f'(p c )), 

{ (f)-\{x-l)/t) if x t <x<\ 

Pi if 1 - y < x < x t 

p r if max(l — y + tv(p c , p r ), 0) < x < 1 — y 

p c if < x < max(l — y + tv(p r , p c ), 0) 

(247) 



p(t,x) 
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for t± < t < r y , where r y = max(t 2 , (y — ^)/ v (p r , Pc)), h is the time at which 
x t — 1, and x t is defined for t > ti by 

x tl = l + tj'fa), x t = v [(fT 1 ((x t - l)/t), Pl ] > v( Pc , Pl ) > (248) 

(j) If p = p c = p*, and f(p{) = f(p r ), then y G [0, 1] is arbitrary, 9 y = +oo, 
T y = max{y/v{pi,p c ), (y - l)/v{p r , p c )), p s (.) = p y (.), 

p c if min(l - y + tv(p h p c ), 1) < x < 1 

pi if 1 - y < x < min(l - y + tv(pi,p c ), 1) 

p r if max(l — y + tf (p c , p r ), 0) < rr < 1 — y 

p c if < x < max(l — y + tv (p r , p c ), 0) 



p(t,x) = < 



(249) 



(k) If f(pi) < /(p r ) and p = p c > p*, then y G [0, 1] is arbitrary ^ 
(1 - y)/v(ip(p r ),<f(pi)), r y is defined below, p s (.) = p h 



Pc 

ip(p r ) 



if < x < max(l — y + tv(p c , ip(pi)), 0) 

if max(l - y + tv(p c , (p(pi)), 0)<x<l-y + tv(ip(pi), <p(p r )) 

p(t,x)={ (fi(Pr) if 1 ~y + tv((p(pi),(p(p r )) < X < 1 -y + tv(<p(p r ),p c ) 

Pc if l-y + tv(ip(p r ),p c ) <x< l + tf(p c ) 

{ (f)-\(x-l)/t) if l + tf(p c )<x<l 

for t < ti :— min(^,t' 1 ), where 

A ■= y/Hv{pr),pc) - f{Pc)\ 

if < x < max(l — y + tv(p c , ip(pi)), 0) 



(250) 



p(t,x) = < 



<f(pi) if 1 -y + tv(p c ,ip(pi)) < x < max(l - y + tv(<p(pi),<p(p r )),0) 

(f{p r ) if 1 - y + tv((p(pi),ip(p r )) < X < X t 

{ (f)-\( x -\)/t) if x t <x<l 

(251) 
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if t[ < t < 9 y , where x t is defined by 

x t[ = l + t[f(p c ), x t = v [(f)-\(x t - l)/t),<p(j>r)] (252) 

with 

v{<p{p r )Mf>l)) < V [V'r\(Xt - l)/t)MPr)] < f'(<p(Pr)) < (253) 

and 

{f)-\{x t -l)/t)e{ip{ Pr ),p c )] 

Pc if < X < tv(pi, (fi(Pr)) 



p(t,x) = < 



<P(Pr) 



if tv(pi,tp(p r )) < X < X t 



(254) 



( ur\(x-\)/t) if Xt <x<i 

if t[ < 6^ < t < t 2 , where t 2 is the time t at which x t = tv(pi, ip(p r )), 



p(t,x) 



Pi 



if < x < x'. 



(255) 



(f')-\(x-l)/t) if x' t <x<\ 



for t[ < Qy < t 2 < t < r y , where x' t is defined by 

< = x t2 , x't = v [UT\« ~ I) ft), Pi] > v{ PhPc ) > (256) 
and t v is the time t at which x' t — 1, 



p(t,x) 



Pi 



if < X < tv(p h (f(p r )) 



(fi(p r ) if tv(pi,<f(p r )) < X < 1 -y + tv(if(p r ),p c 

pc if 1 -y + tv((f(p r ), Pc ) < X < 1 + tf'(p c ) 

[ (f')-\(x-l)/t) if l + tf( Pc )<X<l 



(257) 
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for 9 y < t < min^, t 2 ), where t 2 is the time t at which tv(pi,<p(p r )) 

1 -y + tv{(p{pr),Pe), 



Pi 



if < x < (t - t 2 )v(pi, p c ) + t 2 v(pi, <f(pr)) 



p(t,x) = < 



(258) 



Pc if (t -t 2 )v(pi,p c ) +t 2 v(pi,ip(p r )) < X < 1 +tf'(p c ) 

{ (fTHix-V/t) if 1 + tf(p c ) < x < 1 
if 9 y < t 2 < t < t 2 , where t 2 < t[ is the time t at which 

1 + tf'(pc) = (t- t 2 )v(pi, p c ) + t 2 v(p h (f(pr)) 

pi 



if < x < x'l 



p(t,x) 



(259) 



(f)-\(x-l)/t) if x'l<x<\ 



for t 2 < t < r y , where x" is defined by 

4. 2 = ~L + t' 2 f'(Pc), x" t = v[(f)- 1 ((x';-l)/t),p l ]>v(p l ,p c )> (260) 
and r y is the time t at which x" t — 1, 

P« if < x < tv(p h ip(p r )) 

(p(p r ) if tv(p U ip(p r )) < X < X t 

{ U')-\{x-l)lt) if a;i<a;<l 



(261) 



if 9 y < t\ < t < t 2 , where x t is defined by (I252p . and t 2 is the time t at which 

X t = tv(pi,ip(p r )), 



p(t,x) 



Pi 



if < x < x 



(f')-\(x-l)/t) if xf ] <X<1 



(3) 



(262) 



for t 2 < t < t v , where xf^ is defined by 



,(3) 



x:,/ = x t >>, x^ t = v 



(3), 



UT 1 ^ ~V/t),Pi >r{p,.p,)>() (263) 
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and r y is the time t at which x" — 1. 



(1) If f(pi) < f(p r ) and p = p c < p*, then y G [0, 1] is arbitrary, $ = 
(y - l)/v(ip(p r ), ip(pi)), r y is defined below, p s (.) = p h 

' {fY\xlt) if 0<x<tf'(p c ) 
Pc if tf(p c )<x<l-y + tv(p c ,(p(pi)) 

p(t,x) = < (p(pi) if 1 -y + tv(p c ,ip(pi)) < x < 1 -y + tv(<p(p r ),<p(p l )) 

(p(p r ) if 1 -y + tv(ip(p r ),ip(pi)) < x < min(l,l -y + tv(<p(p r ), p c )) 

k Pc if min(l, 1 - y + tv(<p(p r ),p c )) < x < 1 

(264) 

for t < h := (1 - y)/[f'(p c ) ~ v(p c , <p(pi))}, 
{ {f)-\xlt) if 0<x<^ 

<p(p«) if x t < x < 1 - y + tv((p(p r ),(p(pi)) 

(f{p r ) if 1 -y + tv(<p(p r ),ip(pi)) < x < min(l,l -y + tv(ip(p r ), p c )) 

y Pc if min(l, 1 - y + tv(<p(p r ),p c )) < x < 1 

(265) 

for t 1 < t < t 2 , where t 2 is the time t at which x t — 0, with x t defined by 

= tlf'(Pc), Xt = V[lfi(pi), {f'Y l (x t /t)\ < V(p c , lf(pi)) < V(lf(p r ), lf(pi)) 

(266) 

(p(pi) if xt < x < 1 -y + tv((p(p r ) : (p(pi)) 
p(t,x) = < <p(p r ) if 1 -y + fu(<p(p r ), ip(pi)) < x < min(l, 1 - y + tv(<p(p r ), p c )) 
p c if min(l, 1 - y + tv(ip(p r ), p c )) < x < 1 

(267) 



p(t,x) = < 
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for t 2 < t < 9 y , 

' Pl if o<x<(t-ey)v( Pl ,cp(p r )) 

p(t,x) = < <p(p r ) if 1 -y + tv(ip(p r ),ip(pi)) < x < min(l,l - y + tv(ip(p r ), p c )) 

k p c if min(l, 1 - y + tv(ip(p r ), p c )) < x < 1 



(268) 



if 



y 



< 



v(pi,<p(Pr)) 1 -V(<p(pi),<p(p r )) 1 



v(<p(Pr),Pc)) 1 - V(tp(pi),tp(p r )) 1 

and 9 y <t< r y , with r y = 9 y + (p(p r )), 

' pi if <x < (t-9 y )v(p b (p(p r )) 

p(t,x)=< ip( Pr ) if (t-6 y )v(pi,ip(p r )) < X <l-y + tv(lfi(p r ),p c ) 

^ p c if 1 -y + tv(ip(p r ),p c ) < X < 1 



(269) 



if 



and 



2/ > 



9 y < t < t 3 := (1 - i/) 



V((p(p r ), Pc))' 1 - V^(Pr))" 1 

w(V?(p r , <p(p,)) - ¥»(p r )) 



[u(¥»(/0 r , v(Pi))][^(v(PJ> - V(<p(p r ), C)] 



p(t,x) 



p, if < x < t 3 ^(Pi, y(Pr)) + (* - h)v{pup c ) 

Pc if t 3 v(p h (p(p r )) + (t -t 3 )v(p h p c ) < X < 1 



if 



2/ > 



yjpiivjpr)) 1 -v{y{pi),y{pr)) 1 

v(<p(Pr), Pc))' 1 ~ V(ip(pi), <p(Pr)) _1 

and t 3 <t < r y , where r* is the time t at which 

t 3 v(p h (fi(pr)) + (t- t 3 )v(p h p c ) = 1 



(270) 
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C Explicit minimizers for uniform p(.) when 

Pi > Pr 

Let us assume p(.) = p. Then H195[) yields 

F p q = min[pi,max(p r ,<p(p))], G(0, .) = max[<p(pj), min(<p(p r ), p))\ := p' 

(271) 

C.l The case p r < pi < p* 

In this case, we have 

' p' if x < max(l + t/V),0) 

G(t,a;)=< - if max (! + */' (p')> 0) < x < max(l + tf (<p(pt)), 0) 

k <p(p«) if max(l + tf'(<p(pi)), 0) < x < 1 



Thus 



^(G(t,0 + ))= ^ 



(272) 

f vV) if *<*i := - V/V) 

^[(/O-H-V*)] if *i < * < *2 := -V/'MpO) (273) 
^ Pi if t > t 2 

¥>(G(M~))=P< 



(274) 



From this we obtain: 

First case: ip(pi) < p. Then p > p', and 

p if < x < max(0, 1 + tf'(p)) 



p(t,x) 



for £ < t 2 , 



(/THOe- 1)A) if max(0,l + t/'(p)) < x < 1 



(275) 



p(t,x) 



ft 



if < x < 



(276) 



(f)^(^-l)A) if x t <x<l 
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for t > t 2 , where x t is given by 

x t2 =0, x t = v[p h (fY^ixt - I) ft)] > v[p h (fT^-l/t)} > for t > h 

(277) 

so that Xt reaches 1 in finite time. Note that at time t%, the rarefaction wave 
issued from 1 has value tp(pi) at x = 0, and value < ip(pi) ahead, which is 
why a shock starts propagating from x = towards x — 1. 

Second case: p* < p < ip(pi). Then p' = <p{pi), <f[G(t, + )] = pi for all 
t > 0, and 



p(t,x) 



Pi if < x < tv(pi,p) 

p if tv(p l ,p)<x<l + tf(p) (278) 



[ if)-\(x-l)/t) if l + t/'(p)<x<l 
fort<t 3 := [v(p l ,p)-f(p)}-\ 

{Pi if < x < y t 

(279) 
(f)-\( x -l)/t) if y t <x<l 

where y 4 is given by 

y t3 = t 3 v( Pl ,p), y t = v[pi, {f)-\{y t - l)/t)\ > v{p h p) > (280) 
and thus reaches 1 in finite time. 

Third case: p < p* < <p(pi) < <p(p r )- Then again p' = <p(pi) and (p[G(t, + )] = 

Pi, 

{pi if < x < min[l,tf (pi, p)] 
(281) 
p if min[l, tv(p t , p)] < x < 1 

C.2 The case p r < p* < pi 

First case, p > p* . Then p' > p*, 

(/') _1 ((a: - I) /t) if max[0, 1 + tf'(p')] < x < 1 
G(t,x)={ (282) 
p' if < x < max[0, 1 + tf(p')} 
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and thus 

f <p<jj) if t<t x := -1//V) 

¥>[G(t,0 + )] = { (283) 

I OTH-VO if f >*i 

^[G(t, I")] = (284) 

and 

( {f)-\{x-l)/t) if max[0,l + t/V)]<z<l 
p(t,a;) = { (285) 
[ p' if < x < max[0, 1 + tf'(p')} 

Second case, p < p* . Then p' < p*, 

( (/')~V/0 if O < x < min[M/V)] 
G(t,x) = l (286) 

[ p' if min[l, */'(//)] < x < 1 

and thus 

f y?(p') if * <*i := 1/fV) 

<p[G(t,l-)] = { (287) 

i (/rw*) if *>*i 

0+)]=p* (288) 

and 

(f'r^/t) if 0<x<min[l,i/V)] 
p(t,x)=<( (289) 

p' if min[l,i/'(^')] < rr < 1 

C.3 The case p* < p r < pi 

This case is symmetric to p r < pi < p*, with the roles of the boundaries 
exchanged. Here we have 

' ft if min(M/' (//)) < x < 1 

G(*,x) = < (/'^((a; - l)/t) if min[l,^(p r ))] <x<min[l,*/V)] 
k <p(p r ) if < x < rnin[l,i/'(<p(p r ))] 

(290) 
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Thus 



y?(G(t,0 + )) = p r 
{ <p(pf) if t< tl := 1//V) 



¥>(G(f, 1-)) = I 



V [{f)-\l/t)\ if ti < t < t 2 := l/f(<p(pr)) 



(291) 
(292) 



From this we obtain: 



if t > t 2 



First case: p < (p(p r ). Then p < p', and 

U')-\x/t) if 0<x<mm(l,tf(p)) 



p(t,x) 



for t < t 2 , 



p(t,x) = 



if min(l,£/'(p)) < x < 1 



(f')-\x/t) if 0<x<x t 



(293) 



(294) 



o r if xt < x < 1 

for t > t 2 , where x t is given by 

x t2 =1, x t = v[p r , {f)-\x t /t)\ < v[p r , {f)-\l/t)\ < for t > t 2 (295) 

so that x t reaches in finite time. Note that at time t 2 , the rarefaction wave 
issued from has value <p(p r ) at x = 0, and value > <p(p r ) behind, which is 
why a shock starts propagating from x — 1 towards x — 0. 



Second case: <p(p r ) < p < p* ■ Then p' = (p(p r ), <p[G(t,l~ 
t > 0, and 



p r for all 



p(t,x) = < 



p r if 1 + tV (p r , p) < X < 1 

p if tf'(p) <X< l+tv( Pr ,p) 

( ur l {x/t) if o<x<t/'(p) 

for i <h ■= [f (p) -v(p r ,p)]-\ 

p r if y t < x < 1 

U')-\x/t) if 0<x<y t 
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(296) 



p(t,x) = 



(297) 



where y t is given by 

y t3 = 1 + t 3 v(p r ,p), y t = v\p r , tf')-\Vt/t)\ < v(p r ,p) < (298) 
and thus reaches in finite time. 

Third case: <p(pi) < <p(p r ) < P* < P- Then again p' = <p(p r ) and <p[G(t, l - )] = 

Pr, 

{p r if max[0, 1 + tv(p r , p)] < x < 1 
(299) 
p if < rr < max[0, 1 + tv(p r , p)\ 
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